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Abstract 


Five  assumptions  regarding  a  first  order  model  developed  to  calculate  dose  rate  at  a 
detector  above  a  ground  distributed  fallout  field  are  analyzed.  The  omission  of  scattering  is 
relaxed  by  using  an  integral  method  of  successive  scatters.  The  assumption  of  using  0.70MeV 
as  the  single  average  photon  energy  to  represent  an  entire  fallout  distribution  is  analyzed  using 
the  Oak  Ridge  National  Laboratory  Isotope  Generator  (ORIGEN)  Fallout  Analysis  Tool 
provided  by  Oak  Ridge  National  Laboratory.  A  normalized  Gaussian  distribution  is  used  to 
analyze  the  effects  of  non-uniform  activity  density  on  the  source  normalization  constant 
(SNC).  The  effects  of  time  on  the  SNC  are  also  examined  by  the  ORIGEN  Fallout  Analysis 
Tool  in  order  to  determine  if  a  separate  time  correction  factor  is  needed.  Finally,  previous 
research  accomplished  by  Herte  regarding  self-shielding  from  terrain  roughness  is  summarized. 

Using  a  weighted  average  of  fallout  energy  distribution,  an  average  build  up  factor  of 
1.1  and  an  average  fallout  field  photon  energy  of  0.81MeV  is  calculated.  By  comparing  these 
values  with  a  multi-group  calculation,  they  are  found  to  be  an  accurate  approximation  of  the 
fallout  field  and  result  in  a  25%  increase  to  the  first  order  model  SNC.  The  effects  of  non- 
uniform  activity  density  are  found  to  be  insignificant;  however,  the  effects  of  time  on  the  SNC 
are  found  to  warrant  an  additional  time  correction  factor.  In  summarizing  Herte’s  results, 
terrain  roughness  is  shown  to  produce  between  2-5%  self-shielding  for  soil.  This  value  is 
highly  variable  depending  on  the  terrain  considered,  and  is  not  incorporated  in  the  updated 
SNC  value. 
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SOURCE  NORMALIZATION  CONSTANTS  FOR  GROUND  DISTRIBUTED 


FALLOUT  FIELDS 

I.  Introduction 

Perhaps  one  of  the  most  lingering  attributes  of  a  nuclear  weapon  burst  is  fallout. 
Potentially  covering  vast  areas,  fallout  fields  are  capable  of  producing  many  biological  effects 
ranging  from  acute  radiation  syndrome  to  cancer.  Whether  conducting  a  forensic  analysis 
mission  or  weighing  the  risks  involved  with  rescue  efforts,  the  dose  rate  experienced  while 
traversing  a  fallout  field  is  of  paramount  concern  to  all  involved. 

The  fallout  from  a  nuclear  weapon  burst  is  distributed  across  the  ground  downwind  of 
the  burst.  The  activity  density  of  the  fallout  is  then  converted  to  dose  rate  through  use  of  the 
source  normalization  constant  (SNC).  An  often  used  approximation  of  the  SNC  involves  five 
simplifying  assumptions  that  are  analyzed  at  length  in  this  thesis  with  the  intent  to  improve  the 
approximate  SNC  value.  This  analysis  begins  with  a  reproduction  of  the  first  order  model 
(FOM)  proposed  by  Bridgman  that  is  used  to  calculate  the  current  SNC  value. 

Background 

As  developed  in  Chapter  12  of  The  Introduction  to  the  Physics  of  Nuclear  Weapons  Effects 
(Bridgman,  2001:414),  consider  a  detector  positioned  one  meter  above  a  fallout  field  as 
depicted  in  Figure  1 . 
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Figure  1.  Detector  at  height  h  (lm)  above  a  fallout  field 


The  dose  rate  due  to  ground  distributed  gamma  ray  activity  at  a  given  detector  position 
is  given  by  the  expression: 

(  u  ITzv'h 

(hv)  dArea  d(hv) 

V  P  J  det  (2) 


D=  J  J  Agd[r,t,hv]  BUF[hv] 


hv= 0  Area 


Atts2 


Where, 


D  is  the  dose  rate  in 


J 


kg  -  sec 


at  the  detector ; 


Asd  [rP  ,  hv]  is  the  activity  density 


in 


Bq 


nr  -  J  of  the  y  spectrum 


on  the  ground  at  a 


location  r  measured  from  the  detector; 

BUF[hv]  is  the  energy  dependent  build  up  factor  due  to  scattering  of  gamma  rays  both 
in  the  air  and  ground  [dimensionless]; 

s  is  the  slant  range  from  a  differential  area  (dArea)  to  the  detector  [m] ; 


2 


t  is  time  in  [sec] 


r  =  \ls~  —  1  is  the  ground  distance  from  a  point  just  under  the  detector  to  a  differential 
area  ( dArea )  in  [m] ; 

jUt  [/tv]  is  the  energy  dependent  macroscopic  cross  section  for  attenuation  for  the 
medium  through  which  the  gamma  rays  are  traveling  [/??  '] ; 


p 


is  the  energy  and  material  dependent  total  mass  absorption  coefficient1  for 


det 


the  detector  in 


m 


kg 


hv  is  the  energy  per  disintegration  in  [/] ; 


dArea  is  the  differential  ground  area  [ m  ]  and 


d  ( hv )  is  a  differential  gamma  ray  spectmm  energy  in  [/  of  the  y  spectrum ] 


1  Total  mass  absorption  coefficient  includes  absorption  due  to  Photoelectric  and  Compton  events 
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By  nature,  the  gamma  ray  spectmm  of  a  fallout  field  is  discrete  (See  Figure  2). 


Figure  2.  A  subset  of  the  gamma  ray  spectrum  for  a  fallout  field 
generated  by  the  ORIGEN  Fallout  Analysis  Tool  for  a  lkt  90/10  U- 
235/U-238  fission  at  one  hour. 


This  forces  the  modification  of  equation  (1)  by  replacing  integration  across  all  energies  with  a 
discrete  summation  of  all  gamma  ray  energies  present  at  time  t . 

G(t)  °°  .-A  L*v„Js  (  ,,  r/„,,  i\ 

I. 

g= 1  Area 


u \t)  w 

=  Z  j  (^[^^])  BUF[hvg]- 


-Mdh^s  (  juJhvV 

8  (hv)  dArea 


4  7TS" 


V  P  J  det 


(2) 


Where, 

g  represents  individual  gamma  ray  energies  and  G(t )  the  total  number  of  gamma  ray 

energies  which  varies  with  time; 

hvg  represents  the  energy  of  the  gamma  ray  g  ; 


(Agd[rJ])g  is  the  activity  due  to  gamma  rays  in  group  g  ; 

jUt  and  jUa  are  the  macroscopic  cross  sections  for  total  attenuation  and  total  absorption 


in  ym  1  j  at  the  group  energy  hvg  and 
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BUF  is  the  energy  build  up  factor  for  gamma  rays  at  energy  hvg . 

Simplifying  Assumptions 

Evaluating  equation  (2)  for  all  gamma  ray  energies  present  at  a  given  time  in  a  fallout 

field  is  prohibitively  challenging,  thus  Bridgman  makes  several  simplifying  assumptions: 

1.  He  ignores  scattering  in  the  ground/ air  and  only  calculates  the  direct  contribution  of 
gamma  rays  to  the  detector. 

a.  This  implies  the  build  up  factor  (BUF)  is  one. 

b.  This  also  implies  that  any  gamma  rays  directed  downward  (towards  the 
ground)  never  return  to  the  detector  —  meaning  he  has  immediately 
discounted  half  of  all  gamma  rays  emitted. 

2.  He  replaces  the  entire  discrete  spectrum  of  gamma  rays  with  a  single  average  gamma  ray 
energy  of  0.70MeV  at  a  time  of  one  hour. 

a.  This  eliminates  the  need  to  sum  across  all  energy  groups  as  it  effectively 
reduces  equation  (2)  to  a  one-group  calculation. 

3.  He  evaluates  activity  at  one  hour  and  utilizes  the  Way-Wigner  approximation  to  calculate 
activity  at  all  other  times.  Since  not  all  activity  may  be  deposited  on  the  ground  at  a 
given  time,  Bridgman  invokes  the  concept  of  “unit-time  reference”  for  activity  (and 
inherently  dose  rate).2 

4.  He  assumes  the  activity  density  directly  under  the  detector  at  a  position  (X,Y)  is 
uniform  out  to  infinity.  That  is,  he  assumes  that  activity  density  along  r  does  not  vary, 
but  instead  takes  on  the  value  at  the  foot  of  the  detector. 

a.  This  also  implies  that  he  assumes  an  unfractionated  sample  which  means  the 
gamma  ray  spectmm  does  not  vary  along  r  . 

b.  Furthermore,  the  Way-Wigner  approximation  comes  into  question  for 
subsets  of  the  total  fission  product  set  (volatile /refractory). 

5.  He  assumes  the  ground  is  perfectly  flat  with  no  surface  roughness  to  impede  gamma  rays 
emitted  toward  the  detector. 


2  Reference  The  Effects  of  Nuclear  Weapons  (Glasstone,  427)  for  additional  information  on  unit-time  reference 
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The  First  Order  Model 


Given  an  air  equivalent  detector,  these  assumptions  lay  the  foundation  for  a  first  order 
model  as  depicted  in  the  development  of  the  following  equations 


D[X,Y,l]  =  A[X,Y,l] 


(*v)4 

0 


e-Mhvavg]airS 


V  P  JaU 


4  xs~ 


2  K  r dr 


(3) 


Since,  s2  =  r2  +h  equation  (3)  becomes: 


6[v , r, !] = a[ v, r, i](  [>i>'"8 1 1  (Av)  ij 


co  -fit[hv  ]alrs 


V  p  ) 


2  ■ 


-t/.v 


(4) 


Where  the  integral  is  recognized  as  an  exponential  integral  of  the  first  kind: 


D[X,Y,l]  =  A[X,Y,l] 


Pq[hvavg\ 

V  p  J 


(5) 


Where, 


Z)[X,  F,  1]  is  the  dose  rate  at  a  given  location  ( X ,  F)  at  a  time  of  one  hour  post 


detonation  in 


J 


kg  -  sec 


A[X,  Y,  1]  is  the  activity  density  at  the  given  location  (X,  F)  at  a  time  of  one  hour  in 


Bq 


m 


and 


Ex  is  the  exponential  integral  of  the  first  kind  with  an  argument  of  /ut . 


3  The  equivalent  equations  in  the  2  Aug  2011  revision  of  Introduction  to  the  Physics  of  Nuclear  Weapons  Effects  (Bridgman)  are  12-52a, 
12-52b,  12-52c. 
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premise  is  that  dose  rate  may  be  obtained  through  multiplying  an  activity  density  by  a  constant. 

Research  Objective 

The  objective  of  this  research  is  to  analyze  on  the  quality  of  assumptions  1-5.  Once 
this  assessment  is  made,  any  quantitative  impacts  will  be  combined  to  produce  an  updated 
SNC  for  a  ground  distributed  fallout  field  at  a  time  of  one  hour.  This  analysis  begins  in 
Chapter  II  by  examining  the  impact  of  the  “no  scattering”  assumption. 

4  This  value  is  slightly  different  from  the  value  of  6.33x1  CH  6  developed  by  Bridgman  due  to  small  differences  in  the  values  used 
for  cross  section  values. 
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II.  The  Technique  of  Successive  Scatters 


One  method  for  determining  the  contribution  of  scattering  to  the  detector  is  the 
method  of  successive  scatters  (Bridgman,  2001:209).  In  this  technique,  the  overall  scattering 
contribution  is  broken  up  into  flights  delineated  by  how  many  scatters  occur  between  emission 
and  arrival  at  the  detector.  For  example,  the  term  first  flight  refers  to  gamma  rays  that  reach  the 
detector  without  scattering  (otherwise  known  as  direct  flight).  Gamma  rays  scattered  once 
before  reaching  the  detector  are  referred  to  as  second flight ,  and  so  on  with  each  progressive 
flight.  Using  the  direct  flight  contribution  as  a  baseline,  each  successive  flight’s  contribution  is 
added  until  the  relative  contribution  from  the  n-th  flight  falls  below  an  arbitrary  error  tolerance. 

Depending  on  how  many  flights  are  calculated,  this  technique  can  become 
computationally  intensive.  However,  the  primary  advantage  of  this  integral  method  is  its 
flexibility  with  addressing  assumptions.  More  specifically,  it  continues  to  be  valid  as  other 
assumptions  are  relaxed  over  the  course  of  this  thesis.  This  makes  the  technique  of  successive 
scatters  an  ideal  choice  in  relaxing  the  various  assumptions  associated  with  the  FOM, 
beginning  with  assumption  #1:  the  omission  of  scattering. 

*Note:  Scattering  is  the  only  assumption  relaxed for  the  following  development.  All  other  assumptions  for  the 
FOM  are  in  effect  unless  otherwise  annotated. 

Direct  Flight  Contribution 

Applied  to  the  fallout  field,  we  begin  with  the  direct flight  contribution  to  the  detector 
located  at  (X,  Y )  and  at  a  time  of  one  hour.  Since  this  flight  accounts  for  gamma  rays  that 
reach  the  detector  without  scattering,  its  development  is  identical  to  that  of  the  FOM: 


D[X,Y,l\  =  A[X,Y,l\ 


P 


air 


s 


(4) 


Second  Flight  Contribution  from  Air  Scatter 

To  calculate  the  contribution  of  once  scattered  gamma  rays,  we  begin  by  imagining  a 
differential  volume  located  at  some  other  location  {X’ ,  Y')  at  a  time  of  one  hour  and  at  a 
height  h'  (see  Figure  3). 


Figure  3.  Representation  of  a  differential  volume  as  a  scattering  source. 


This  differential  volume  receives  a  direct  flight  contribution  from  the  ground 
distributed  fallout  field.  In  this  respect,  the  propagation  of  gamma  rays  to  this  differential 
volume  is  identical  to  the  previous  direct  flight  development  with  one  exception:  instead  of 
being  absorbed  (as  was  the  case  for  the  detector),  the  gamma  rays  are  scattered  both 
coherendy  (Rayleigh)  and  incoherently  (Compton)  depending  on  cross  sections.  For  now,  we 
will  assume  the  scatter  is  isotropic  and  incoherent.  This  concept  is  illustrated  by  the  equation: 


(7) 
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*  Not e:  Dependencies  have  been  removed  to  simplify  presentation  of  the  equations 

Where, 

hvo  is  the  pre-scattered  gamma  ray  energy  in  [,/] ; 
ju°  is  the  scattering  cross  section  at  hvo  for  air  in  [m  1  ] ; 

//“  is  the  attenuation  cross  section  at  hvo  for  air  in  [m  '] ; 
h'  is  the  height  of  the  differential  volume  in  [ml ; 

s'  is  the  slant  range  from  the  differential  volume  to  a  differential  area  on  the  ground  in 
[ml  and 


S  is  volumetric  scattering  source  within  the  differential  volume  in 


Bq 


m 


As  viewed  from  the  detector,  this  differential  volume  is  identical  to  a  volumetrically 
distributed  source  in  air  that  is  emitting  at  an  activity  S  .  Furthermore,  this  is  only  one  of  an 
infinite  number  of  differential  volumes  in  the  space  around  the  detector.  As  such,  the 
contribution  of  each  differential  volume  must  be  calculated  in  the  form  of  a  spatial  integral: 


Dmr  = 

1J SF 


p , 


CO 

(Av>  j  S, 


-WairS" 


4  jtr 


n  2 


-dVol 


(8) 


Where, 


Dfp  is  the  second  flight  contribution  to  dose  rate  in 


J 


kg  -  sec 


hv'  is  post-scatter  gamma  ray  energy  in  [/] ; 
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— 1 -  is  the  total  mass  absorption  coefficient  of  the  detector  for  gamma  rays  at 

p  J  . 

'  /  air 


energy  hv'  in  -  ;  *Note:  here  an  air  equivalent point  detector  is  assumed 

kg 


p\  is  the  attenuation  cross  section  for  air  at  hv'  in  1 J  ; 

r"  is  the  ground  range  from  directly  beneath  the  detector  to  directly  beneath  the 

differential  volume  in  [m] ; 

5ffis  the  slant  range  from  the  detector  to  the  differential  volume  in  [/n]  and 
dVol  is  the  differential  volume  defined  as  2 nr"dr"dh'  in  j  . 

Similar  to  the  FOM  development,  a  change  of  variables  must  be  accomplished  to 
homogenize  the  variable  of  integration.  While  temporarily  holding  h'  constant: 


s2  =  r2  + 1 h  -  h!\  and  r"dr"  =  s"ds" 


Also  associated  with  the  change  of  variables  is  a  respective  change  in  limits  of  integration: 


r"  =  0  s"  =  \h  —  h' 


When,  1  „  „ 

r  =oo  s  =  oo 


After  executing  the  change  of  variables  and  once  again  allowing  //'  to  vary,  the  result  is: 


A  air  _  Pa 


1  ”  ”  p-WW 

|  - -i-ds'dh' 


n  7J  J  v 

r  J  air  ^  0  \h-h’\  A 


The  absolute  value  in  the  limit  of  integration  necessitates  further  manipulation  by 
breaking  up  the  integral  into  two  parts  based  on  the  relationship  between  h  and  h' . 


Physically,  the  integral  is  broken  into  the  regions  illustrated  in  Figure  4: 
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Figure  4.  Regions  II  and  III  representing  the  space  below  and 
above  the  detector,  respectively. 


(10) 


(11) 


After  the  application  of  another  change  of  variables  for  each  region,  the  total  contribution  to 
dose  rate  from  once  scattered  gamma  rays  in  air  is: 

(12) 


12 


Second  Flight  Contribution  from  Ground  Scatter 

Given  an  isotropic  emission  of  gamma  rays  from  the  ground  distributed  source,  air 
scatter  only  addresses  half  of  the  emissions.  The  other  half  (which  are  ignored  in  the  FOM) 
are  gamma  rays  emitted  directly  into  the  ground  which  are  absorbed,  scattered  into  the  ground, 
or  scattered  back  into  the  air.  Figure  5  illustrates  the  contribution  of  fallout  to  a  differential 
volume  in  ground.  Again,  all  scatters  are  still  assumed  to  be  isotropic  at  this  point. 


Figure  5.  A  differential  volume  representing  a  scattering  source  in  the 
ground. 


The  contribution  made  to  the  detector  from  gamma  rays  once  scattered  in  the  ground 
begins  similar  to  the  air  scatter  development  in  equation  (7)  with  the  modeling  of  a  source 
term  describing  the  propagation  of  gamma  rays  to  a  differential  volume: 


sgd=*gd 


(13) 
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The  only  differences  from  the  previous  source  term  development  are  the  scattering  and 
attenuation  cross  sections  now  reflect  the  change  in  medium  from  air  to  ground. 

Also  similar  to  the  air  scatter  example  is  the  development  of  the  dose  rate  expression. 
Once  scattered  from  the  differential  volume,  the  journey  of  gamma  rays  to  the  detector  must 
be  broken  into  two  segments:  one  for  ground,  and  one  for  air.  This  is  illustrated  in  Figure  6: 


Figure  6.  A  gamma  ray  scattered  within  the  ground  must  travel  through 
both  ground  and  air  as  it  proceeds  to  the  detector. 


bgd  = 

^ SF 


p , 


(*vyl 


J  V 


-(Mi)gd  *2'-  («’)«,>  *1' 


-ds" 


h+h' 


Where, 


si"  is  the  slant  range  distance  covered  in  the  air  in  [m]  and 
s 2"  is  the  slant  range  distance  covered  in  the  ground  in  [/«] . 

After  applying  some  trigonometry  and  a  change  of  variables,  the  end  result  is: 


(14) 
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D«=c]ElUM;)'l(u-h) 


E, 


M, 


gd 


(u-h)  +  [/u'\  h 

\  /  air 


\du 


(15) 


Finding  Post-Scatter  Energy 

Before  a  build  up  factor  (BUF)  is  calculated,  a  few  hurdles  remain  in  obtaining  a  value 
for  hv'  contained  in  the  constant  term  (C)  of  equations  (12)  and  (16).  Using  the  three  basic 
conservation  equations,  The  A  tomic  Nucleus  gives  us  the  following  relationship  for  Compton 
scatter  (Evans,  1955:675): 

hv’  _  1 

hv0  l  +  a(l-cos<9) 


Where, 


hv'  is  the  post-scatter  energy  in  [MeU ] ; 

hv0  is  the  incident  energy  in  [MeU ] ; 
hv 

a  is  the  value  - V  and 

m0c 

1 9  is  the  photon  angle  of  deflection. 

This  relationship  illustrates  the  direct  correlation  between  angle  of  deflection  and  post¬ 
scatter  energy.  Accounting  for  all  possible  angles  and  post-scatter  energies  for  all  incident 
energies  is  extremely  challenging,  thus  a  simplifying  approximation  is  made  in  using  an  average 
angle  of  deflection  to  determine  an  average  post-scatter  energy.  This  is  accomplished  using  a 
derivation  of  the  Klein-Nishina  formula  to  find  the  angular  distribution  of  Compton  scattered 
photons  (Evans,  1955:690): 


d  Qo-) 

dS 


^2;rsin,, 

dQ. 


(17) 
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Using  an  expanded  version  of  equation  (17),  a  plot  is  generated  for  the  0.70MeV  case  (shown 
in  Figure  7): 


Figure  7.  Plot  representing  the  number  of  scattered 
photons  per  unit  angle 


By  integrating  equation  (1 7)  from  0  to  n,  the  entire  area  of  the  lobe  displayed  in  Figure 
7  is  calculated.  Using  this  value  as  a  normalization  factor,  a  weighted  average  of  the  number  of 
photons  per  angle  is  taken  to  obtain  an  average  angle  of  scatter.  This  average  angle  is  used  in 
equation  (16)  to  obtain  an  average  energy  ratio,  which  is  then  multiplied  by  the  incident  energy 
to  get  an  average  post-scatter  energy. 

Ground  Non-isotropic  Adjustment  Term 

Even  though  Figure  7  is  evidence  of  the  contrary,  up  to  this  point  it  is  assumed  that  all 
scatters  are  isotropic.  As  shown  in  Figure  8,  this  assumption  clearly  breaks  down  as  energy 
increases  and  gamma  rays  are  preferentially  scattered  in  the  forward  direction. 
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Figure  8.  Klein-Nishina  plot  from  The  Atomic  Nucleus  (Evans,  1955:683) 


Accounting  for  all  possible  scattering  angles  using  an  integral  method  is  challenging; 
however,  a  rough  approximation  can  be  made  in  the  case  of  ground  scattered  gamma  rays  to 
further  refine  the  model.  Given  a  fallout  field  distributed  over  a  smooth  surface,  the  range  of 
possible  scattering  angles  a  gamma  ray  could  achieve  in  ground  near  the  detector  in  order  to 
scatter  back  into  the  air  is  90°-270°.  Therefore  at  any  given  energy,  the  ratio  of  the  backscatter 
region  area  (90°-270° )  to  the  area  of  the  entire  scattering  range  (0°-360°)  gives  a  rough  estimate 
of  the  fractional  portion  of  ground  scatter  that  returns  to  the  air. 

The  obvious  flaw  in  this  adjustment  term  is  that  gamma  ray  emissions  from  large 
horizontal  distances  may  graze  the  ground,  and  forward  scatter  towards  the  detector.  This  is 
why  the  stipulation  “ near  the  detected’  is  necessary  when  describing  the  development  of  this 
term.  This  value  is  designated  as  the  Ground  Non-isotropic  Adjustment  Term  (GNAT). 
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This  temporarily  assumes  that  scatter  in  ground  is  completely  isotropic  —  but  for  purposes  of 
this  additional  refinement  to  the  BUF,  the  approximation  is  sufficient.  Mathematical 
development  of  the  contributions  of  this  additional  term  is  identical  to  that  of  the  second  flight 
air  contribution  derivation  from  a  previous  section  (see  page  9). 
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III.  Scattering  Contributions  to  the  Build  Up  Factor 


The  Modeling  Environment 

To  obtain  the  results  discussed  in  this  section,  various  assumptions  about  air  and 
ground  composition  are  made.  The  reader  is  cautioned  that  every  number  reported  in  this 
section  is  not  absolute,  and  must  be  associated  inseparably  with  these  assumptions  and 
approximations.  In  other  words,  consider  all  values  to  have  qualitative  error  bounds 
associated  with  the  various  assumptions  used  in  this  model. 

In  all  models,  ground  is  assumed  to  be  homogeneous  in  both  composition  and  density. 
As  ground  composition  varies  significantly,  a  simplification  is  made  by  assuming  a  chemical 
composition  of  100%  Si02.  Also  varying  with  composition,  location,  compaction,  and  several 
other  factors  is  soil  density.  Density  values  range  anywhere  from  1000-2000  kg/m3  for  sandy 
soils,  so  a  simplified  value  of  1922  kg/ m3  is  arbitrarily  assumed.  Finally,  the  ground  is  assumed 
to  be  a  perfectly  smooth  and  level  plane  (assumption  #5  of  the  FOM). 

Applied  to  all  models,  air  is  assumed  to  contain  78%  N2  /  21%  02  /  1%  Ar. 

Following  the  International  Standard  Atmosphere,  a  temperature  and  density  of  1 5°C  and  1 .23 
kg/ m3  for  sea  level  air  are  observed.  These  parameters  are  taken  as  universal  and  homogenous 
at  all  locations  and  altitudes. 

Fallout  distribution  is  assumed  to  be  homogenous  and  unfractionated  in  all  cases.  As 
stated  in  the  FOM,  the  activity  density  directly  under  the  detector  is  assumed  to  be  uniform  at 
all  locations.  The  variation  of  activity  density  over  the  distribution  of  fallout  is  analyzed  in 
Chapter  V. 

For  cases  involving  a  single  average  gamma  ray  energy,  the  source  is  assumed  to  emit 
at  the  specified  average  energy  in  an  isotropic  fashion.  In  evaluation  of  the  entire  fallout 
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energy  spectrum,  activity  weighted  values  are  used  to  describe  the  fractional  energy 
contribution  of  each  individual  energy  group.  These  fractional  contributions  are  assumed  to 
be  constant  at  all  locations. 

As  previously  stated,  all  scattering  events  are  assumed  to  be  isotropic  unless  otherwise 
noted.  This  assumption  is  obviously  incorrect  as  shown  in  the  Klein-Nishina  plots  (Figure  8); 
however,  it  is  assumed  that  the  overestimation  of  scattering  contributions  for  trajectories 
pointed  away  from  the  detector  will  partially  balance  with  the  underestimation  of  scattering 
events  that  occur  for  trajectories  pointed  at  the  detector.  The  quality  of  this  assumption  is  not 
addressed  and  it  is  listed  as  a  simplifying  assumption. 

Fundamental  to  modeling  all  photon  interactions  with  matter  are  cross  sections.  All 
cross  section  values  were  obtained  from  the  National  Institute  of  Standards  and  Technology 
(NIST)  XCOM  Photon  Database.  For  energies  not  listed  on  the  table,  a  simple  linear 
interpolation  scheme  (on  a  log-log  plot)  is  used. 

Finally,  sample  calculations  for  each  of  the  numerical  values  reported  in  this  chapter 
may  be  found  in  Appendix  A  referenced  by  superscripts  following  each  reported  value  in  this 
chapter. 

The  Case  for  a  Single  Incident  Average  Energy 

As  proposed  in  The  Effects  of  Nuclear  Weapons  (Glasstone,  1977:454),  the  average  energy 
of  all  gamma  rays  at  a  time  of  one  hour  after  detonation  is  0.70MeV.  Using  the  techniques 
previously  discussed  in  Chapter  II,  0.70MeV  gamma  rays  Compton  scatter  down  to  an  average 
energy  of  approximately  0.38MeV’.  At  this  post-scatter  energy,  isotropic  scattered  gamma  rays 


5  See  Appendix  A  (1) 
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in  air  contribute  7.5%6  of  the  dose  rate  as  seen  by  the  detector  -  with  83%7  of  the  air  scattering 
contribution  coming  from  Region  III  (the  region  above  the  detector).  In  addition,  isotropic 
gamma  ray  scatter  from  the  ground  contributes  6.1%8  of  the  total  dose  rate. 

After  making  further  adjustments  for  non-isotropic  ground  scatter  using  the  Ground 
Non-isotropic  Adjustment  Term  (GNAT)  and  adding  the  contributions  from  the  limited 
second  scatter  approximation,  the  final  BUF  value  for  single  average  incident  energy  of 
0.70MeV  is  calculated  to  be  1.1. 9  These  values  result  in  a  source  normalization  constant  of10: 


SNC  for  0.70MeV  and  BUF  of  1 . 1 


7.5xl(T16 


J-m2 
kg  -  sec-  Bq 


This  new  SNC  replaces  the  constant  previously  mentioned  in  equation  (6).  These  results 
indicate  that  scattering  contributes  at  least  10%  to  the  SNC. 

The  impact  of  the  BUF  is  not  constrained  to  a  single  energy,  however.  As  such,  the 
BUF  value  associated  with  a  single  average  gamma  ray  energy  of  0.70MeV  must  be  contrasted 
with  BUF  values  for  the  actual  distribution  of  fallout  energy  to  further  gauge  its  quality  as  an 
assumption. 

The  Case  for  Multiple  Incident  Energies 

As  illustrated  in  Figure  9,  the  previously  discussed  average  gamma  ray  energy  0.70MeV 
falls  within  the  Compton  scattering/ absorption  region  for  air.  In  the  energy  range  of  0.30MeV 
-  l.OMeV,  the  lines  for  total  absorption,  total  attenuation,  and  Compton  scattering  are  fairly 


6  See  Appendix  A  (2) 

7  See  Appendix  A  (3) 

8  See  Appendix  A  (2) 

9  See  Appendix  A  (4) 

10  This  value  is  computed  by  multiplying  the  previous  SNC  of  6.8xl0‘16  by  the  BUF  of  1.1 
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constant  and  remain  in  the  same  order  in  terms  of  the  mechanisms  of  matter  interaction. 
However,  the  region  to  the  far  left  (roughly  0.02MeV  -  O.lOMeV)  relative  cross  sections  have 
steeper  slopes.  The  question  then  becomes,  what  impact  might  this  have  on  the  BUF?  And, 
will  this  impact  be  significant  enough  to  further  change  the  source  normalization  constant? 


Mev 


Figure  9.  Mass  coefficients  for  air  from  The  Atomic  Nucleus  (Evans,  1955:713). 


Utilizing  the  identical  model  and  method  used  for  the  single  average  incident  energy 
case,  multiple  incident  energy  groups  were  analyzed  to  determine  a  BUF  across  the  spectmm 
of  interest.  The  specific  breakdown  of  these  energy  groups  and  the  rationale  for  the  energy 
range  depicted  are  addressed  in  Chapter  IV.  These  results  are  displayed  in  Figure  10: 
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Figure  10.  The  BUF  as  a  function  of  incident  gamma  energies 


As  shown  in  Figure  10,  the  BUF  rises  to  nearly  1.9,  rapidly  declines,  and  then 
asymptotically  approaches  a  value  close  to  1 .0.  This  relationship  is  explained  by  understanding 
the  inherent  limitations  to  the  BUF  value  and  a  brief  analysis  of  cross  sections. 

The  first  observation  is  the  apparent  asymptotic  approach  to  a  BUF  of  1.0  as  energy 
increases.  This  can  be  attributed  to  a  combination  of  cross  section  behavior  and  the  energy 
lost  due  to  a  scattering  event  at  high  energies.  Shown  in  Figure  1 1,  as  energy  increases  the 
amount  of  energy  lost  due  to  a  scattering  event  also  increases. 


11  This  chart  is  the  result  of  the  BUF  development  from  the  previous  section  applied  to  multiple  energies. 
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Average  Scatter  vs.  Incident  Energy 


1.000 


Incident  Gamma  Energy  (MeV) 


Figure  11.  Average  scattered  energy  as  a  function  of  incident  energy12 


The  larger  the  amount  of  energy  loss,  the  greater  the  shift  to  the  left  on  Figure  9  for 
post-scatter  energy  -  which  means  a  greater  change  in  relative  cross  section  values  from  pre  to 
post  collision.  In  the  region  of  interest  for  nuclear  fallout  (O.OMeV  —  3.5MeV),  a  shift  to  the 
left  universally  equates  to  an  increase  in  total  attenuation  cross  section  and  a  decrease  in  mean  free 
path.  This  creates  a  situation  where  high  energy  gamma  rays  have  the  energy  to  travel  far  away 
from  the  detector  before  interacting  with  matter;  however,  once  they  interact  and  lose  a 
significant  fraction  of  their  energy,  their  chances  of  making  the  same  trip  back  to  the  detector 
without  attenuating  are  severely  reduced.  As  such,  the  BUF  is  expected  to  asymptotically 
approach  a  value  close  to  1.0  as  energy  increases. 


12  This  chart  is  an  application  of  Appendix  A  (1)  across  multiple  energies. 
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Figure  1 0  also  exhibits  a  peak  in  BUF  value  associated  with  lower  energies.  Again,  this 
feature  is  explained  by  analyzing  the  cross  sections  featured  in  Figure  9.  In  this  case,  the 
energy  lost  in  low  energy  collisions  is  less  than  1%.  Nonetheless,  even  the  smallest  decrease  in 
energy  results  in  a  relatively  large  increase  in  absorption  cross  section,  yet  only  small  relative 
change  in  scattering  cross  section.  In  this  scenario,  low  energy  gamma  rays  will  tend  to  scatter 
(losing  very  little  energy  in  the  process)  until  absorbed  by  the  detector. 

The  BUF  value  for  the  range  below  0.02MeV  is  associated  with  the  physical  geometry 
of  the  problem  and  is  predicted  to  eventually  decline  to  a  value  of  1.0.  This  concept  is  best 
described  using  Figure  12: 


Figure  12.  Geometric  limitations  of  low  energy  scatters 


Consider  two  identical  gamma  rays  (A  and  B)  emitted  from  the  ground.  If  the 
problem  is  constrained  such  that  a  gamma  ray  must  scatter  once  before  reaching  the  detector, 
there  is  a  limited  volume  of  interaction  (shown  above  as  the  sphere  surrounding  the  detector) 
in  which  a  gamma  ray  can  scatter  and  still  have  a  significant  probability  of  arriving  at  the 
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detector.  In  Case  B,  the  gamma  ray  scatters  outside  this  limited  volume  and  does  not  possess 
the  energy  to  have  a  significant  probability  of  reaching  the  detector.  As  energy  decreases,  the 
proximity  of  the  scattering  event  to  the  detector  must  also  decrease  if  the  gamma  ray  is  to 
reach  the  detector.  Likewise  as  energy  increases,  the  proximity  of  the  scattering  event  may 
occur  at  larger  distances.  The  concept  of  this  limited  volume  of  interaction  in  the  form  of  a 
sphere  is  also  vital  in  describing  the  relative  contributions  of  air  scatter  and  ground  scatter  to 
the  detector  (Figure  13). 


Relative  Scattering  Contributions  to  Dose  Rate  (Isotropic  Scatter) 
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Figure  13.  Relative  contributions  of  scattered  energy  from  air  and  ground13 


-Air 

Ground 


As  shown  in  Figure  13,  air  is  the  dominant  source  in  scattering  contributions  to  the 
detector.  Conceptually,  this  is  understood  by  comparing  the  two  mediums.  In  air,  attenuation 


13  This  chart  is  an  extension  of  the  values  calculated  in  Appendix  A  (4)  across  multiple  energies. 
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cross  sections  are  fairly  low.  This  allows  for  mean  free  paths  of  18  meters14  for  a  0.02MeV 
gamma  ray  and  110  meters  for  a  0.70MeV  gamma  ray.  Mean  free  paths  on  the  order  of  100 
meters  mean  gamma  rays  have  the  ability  to  traverse  much  larger  distances  and  still  reach  the 
detector.  Conversely,  ground  attenuation  cross  sections  allow  for  mean  free  paths  of 
2. lx  10  3  meters  for  0.02MeV  gamma  rays  and  7.8x10  meters  for  0.70MeV  gamma  rays. 
These  values  suggest  that  unless  the  scattering  event  occurs  extremely  close  to  the  surface  - 
the  gamma  ray  will  likely  never  make  it  out  of  the  ground. 

A  less  obvious  observation  is  the  energy  at  which  relative  contributions  from  these  two 
mediums  peak  are  different.  This  is  explained  by  referencing  Figures  11  and  12.  The  air 
contribution  peak  occurs  at  roughly  0.02MeV.  At  this  energy,  the  limited  volume  of 
interaction  described  in  Figure  12  is  in  optimum  balance.  An  increase  in  energy  from  this  point 
of  balance  would  increase  the  limited  volume  of  interaction,  but  that  increase  in  volume  would 
actually  result  in  a  decrease  in  first  scatter  contributions  to  the  detector.  This  decrease  is  because 
of  spherical  divergence  and  the  previously  described  case  of  a  gamma  ray  having  enough 
energy  to  travel  a  large  distance  before  scattering,  but  not  having  enough  energy  to  return. 
Likewise,  a  decrease  in  energy  from  the  point  of  balance  would  result  in  a  decrease  in  the 
limited  volume  of  interaction.  This  smaller  volume  equates  to  a  smaller  quantity  of  gamma 
rays  being  able  to  scatter  and  still  make  it  to  the  detector.  In  essence,  the  air  contribution  peak 
in  Figure  13  is  a  result  of  the  Goldilocks  adage  applied  to  the  volume  of  interaction  -  not  too 
big,  not  too  small,  but  just  right. 

As  it  turns  out,  this  ideal  limited  volume  of  interaction  for  air  is  far  from  ideal  for 
ground.  Although  ground  contributions  are  on  the  rise  at  0.02MeV,  they  are  not  at  a  peak. 


14  Appendix  A  (5) 
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This  is  because  the  aforementioned  volume  of  interaction  only  barely  intersects  the  ground, 
meaning  the  volume  of  air  available  for  interaction  dwarfs  the  volume  of  ground  available  (see 
Figure  14). 


Because  of  the  minor  amount  of  ground  contained  within  the  limited  volume  at 
0.02MeV,  it  is  understandable  that  scattering  contributions  from  the  surface  will  be  relatively 
minor  -  especially  considering  ground  interactions  are  more  dependent  upon  surface  area  than 
volume  as  described  before  with  mean  free  path  comparisons.  Furthermore,  any  ground 
interaction  will  also  need  to  traverse  a  volume  of  air  before  reaching  the  detector.  After  energy 
increases  to  approximately  O.lOMeV,  the  volume  of  limited  interaction  grows  into  the 
previously  mentioned  balanced  condition  for  the  same  reasons  -  only  this  time  conditions  are 
optimized  for  ground  scatter  contributions. 

Another  observation  is  that  ground  scatter  contributions  never  exceed  the 
contributions  of  air  scatter,  but  instead  approach  a  state  of  almost  identical  contribution  at 
higher  energies.  Again,  using  the  concept  illustrated  in  Figure  14  this  makes  sense.  The 
volume  of  air  within  the  limited  volume  of  interaction  will  always  exceed  the  corresponding 
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volume  of  ground  because  the  reference  volume  is  centered  on  the  detector.  This  disparity  in 
volume  that  dominates  low  energies  becomes  less  apparent  as  the  limited  volume  of 
interaction  grows  with  increasing  energy.  As  the  disparity  in  volume  becomes  less 
pronounced,  the  difference  in  relative  air  and  ground  scatter  contribution  also  becomes  less 
pronounced. 

Finally,  a  caution  must  be  observed.  Although  the  mechanics  of  these  interactions 
occur  as  previously  described,  definitively  associating  these  peaks  with  a  specific  gamma  ray 
energy  is  somewhat  misleading.  The  fact  that  the  BUF  peak  occurs  at  roughly  0.02MeV  is 
merely  a  product  of  how  the  information  is  presented  (in  terms  of  incident  energy)  and  how 
many  scatters  are  allowed  (one).  If  the  BUF  is  plotted  as  a  function  of  incident  energy  while 
only  considering  two  scatters  —  the  relative  BUF  peak  would  appear  at  slightly  higher  incident 
energy.  The  effect  of  considering  multiple  scatters  is  predicted  to  broaden  the  peak  observed 
in  Figure  10,  effectively  expanding  higher  BUF  values  across  a  larger  spectmm  of  energy.  This 
fact  further  emphasizes  the  need  to  expand  research  of  this  problem  into  multiple  scatters. 

Regardless  of  how  many  scatters  are  considered,  Figure  10  demonstrates  a  definitive 
energy  dependence  of  the  BUF  term.  Because  of  this  dependence,  the  BUF  value  of  1.1  for  an 
average  energy  of  0.70MeV  may  not  be  representative  of  the  actual  fallout  energy  distribution. 
Depending  on  which  energies  dominate  the  spectmm,  the  higher  BUF  values  associated  with 
lower  energies  may  further  influence  the  source  normalization  constant.  This  issue  is 
examined  in  the  next  chapter. 
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IV.  Challenging  the  0.7MeV  Average  Energy  of  Fallout 


As  previously  discussed,  assumption  #2  in  the  development  of  the  FOM  reduces  the 
full  spectrum  of  energy  emitted  from  the  fallout  field  to  a  single  average  gamma  ray  energy  of 
0.70MeV.  In  Chapter  III,  results  from  Figure  10  show  a  definitive  energy  dependence  of  the 
BUF  value.  Depending  on  the  energy  distribution  of  fallout,  the  higher  BUF  values  associated 
with  other  energies  could  impact  the  calculated  value  for  a  source  normalization  constant  at 
one  hour  post  detonation.  As  such,  the  objective  of  this  chapter  is  twofold:  First,  to  verify  the 
accuracy  of  using  0.70MeV  as  the  average  photon  energy  emitted  from  a  fallout  field  at  a  time 
of  one  hour.  Once  the  average  is  verified,  the  second  objective  is  to  determine  if  using  average 
photon  energy  is  a  reasonable  approximation  in  calculating  a  BUF  and  the  resulting  source 
normalization  constant. 

Using  the  ORIGEN  Fallout  Analysis  Tool  provided  by  Oak  Ridge  National 
Laboratory,  fallout  energy  distributions  for  two  different  fission  types  were  analyzed.  Fissile 
material  compositions  are  taken  as  90%/ 10%  Uranium-235/Uranium-238  respectively,  and 
100%  Plutonium-239.  Both  fuels  are  assigned  a  yield  of  1  Kiloton  and  contain  1.0kg  of 
material11.  The  specific  procedure  and  model  values  used  in  these  simulations  are  addressed  in 
Appendix  B  while  actual  code  output  may  be  found  in  Appendix  C. 


15  1.0kg  of  fissile  material  cannot  produce  a  lkt  yield;  however,  proportions  of  the  isotope  inventory  are  unaffected  by  fuel 
mass.  As  such,  assuming  a  1.0kg  mass  does  not  affect  the  values  obtained  in  this  chapter. 
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Average  Energy  and  Total  Photon  Activity  for  100%  Pu-239  Fission 

The  energy  spectrum  from  O.OMeV  to  5.0MeV  is  broken  up  into  50  energy  groups  of 
equal  energy  (Figure  15).  This  range  accounts  for  —  1 00%  of  all  energy  of  the  distribution  at 
one  hour16.  From  there,  an  activity  weighted  average  is  used  to  calculate  the  average  photon 
energy  and  total  activity  within  the  ORIGEN  Fallout  Analysis  Tool. 


Figure  15.  The  energy  spectrum  of  fallout  from  100%  Pu-239  fission  at  one  hour. 


The  average  photon  energy  for  Pu-239  fission  is  found  to  be  0.79MeV17,  while  the  total  activity 


is  shown  to  be  595  Photon  Megacuries  .  These  results  are  shown  in  Figure  16. 


16  See  the  last  lines  of  Figure  16  to  obtain  the  “not  counted”  portion  of  energy 

17  See  Appendix  B  for  model  parameters 


18  This  value  is  for  gamma  rays  and  x-rays 
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Photon  Source  Definition  -  generated  from  ORIGEN  F71  file 
Filename:  today. f 71 

case  106.  3.600E03  seconds  0.000E00 

total  strength:  2.2051E19  gammas/second 
average  energy:  0.7865306796307987 

discrete  lines:  0.0000E00  gammas/second 
0.00%  of  energy 
average  line  energy:  0.0 

multigroup  bins:  2.2051E19  gammas/second 
100.00%  of  energy 

average  group  energy:  0.7865806796307987 

not  counted:  1.4098E11  gamma3/second 

0.00%  of  the  original  energy 


Figure  16.  Results  from  ORIGEN  Fallout  Analysis  Tool  for  average  energy 
and  total  photon  activity  of  fallout  from  Pu-239  fission  at  one  hour. 


Since  these  results  deviate  from  the  widely  accepted  values  of  0.70MeV  and  530 
Gamma  Megacuries19  in  The  Effects  of  Nuclear  Weapons  (Glasstone,  1977:454),  another 
composition  is  analyzed  to  ensure  that  fuel  composition  is  not  a  dominating  factor  in  the 
values  obtained  for  average  energy  and  total  photon  activity. 

Average  Energy  and  Total  Activity  for  90/10  U-235/U-238  Fission 

The  energy  spectmm  from  O.OMeV  to  3.5MeV  is  broken  up  into  35  energy  groups  of 
equal  energy.  This  range  accounts  for  99.6%  of  all  energy  of  the  distribution  at  one  hour2". 
From  there,  an  activity  weighted  average  is  used  to  calculate  the  average  gamma  ray  energy  and 
total  activity  within  the  ORIGEN  Fallout  Analysis  Tool.  The  average  gamma  ray  energy  is 
found  to  be  0.81MeV21.  This  is  in  contrast  with  the  0.70MeV  average  listed  in  ENW. 

19  Although  this  value  includes  only  gamma  rays,  the  lack  of  X-ray  inclusion  does  not  account  for  the  65  Megacurie  difference. 

20  See  the  last  lines  of  Figure  17  to  obtain  the  “not  counted”  portion  of  energy 

21  See  Appendix  B  for  model  parameters 
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Likewise,  total  gamma  ray  activity  is  shown  to  be  583  Gamma  Megacuries  -  compared  to  the 


530  Gamma  Megacuries  listed  in  ENW.  These  results  are  shown  in  Figure  17. 


Photon  Source  Definition  -  generated  from  0RIGE1J  F71  file 
Filename:  today. f71 

case  106.  3.600E03  seconds  0.000E00 

total  strength:  2.1569E19  gammas /second 
average  energy:  0.8123485155581431 

discrete  lines:  0.0000E00  gammas/second 


0.00%  of  energy 
average  line  energy:  0.0 

multigroup  bins:  2.1569E19  gammas/second 
100.00%  of  energy 

average  group  energy:  0.8123485155581431 

not  counted:  1.9829E16  gammas/second 

0.40%  of  the  original  energy 


Figure  17.  ORIGEN  Fallout  Analysis  Tool  output  for  average  energy  and  total 
activity  at  1  hour  for  a  90/10  U-235/U-238  fission. 


To  convert  from  “total  strength”  to  Gamma  Megacuries,  the  following  relationship  is  used: 


Total  Activity  = 


(19) 


The  differences  in  these  values  from  those  used  within  ENW  ate  of  unknown  origin, 


but  most  likely  can  be  attributed  to  the  use  of  modern  data  sources  (such  as  improved  cross 


sections  and  fission  yield  values)  contained  within  the  ORIGEN  Fallout  Analysis  Tool. 
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Inclusion  of  X-rays  into  the  Model 


In  regards  to  the  source  normalization  constant  (SNC),  gamma  rays  are  not  the  only 
source  of  photons  contributing  to  dose  rate.  These  X-rays  originate  from  capture  and  de¬ 
excitation  of  electrons  by  fission  products.  Another  potential  excitation  mechanism  is  ejection 
of  an  inner  shell  electron.  In  this  case,  the  remaining  electrons  shift  to  fill  the  inner  shell 
vacancy  and  emit  X-rays  in  the  process.  Regardless  of  the  excitation  mechanism,  excited 
electrons  shifting  to  a  ground  state  in  elements  containing  more  than  six  protons  cause  the 
emission  of  X-rays. 

X-rays  also  contribute  to  the  dose  rate  seen  by  the  detector.  As  such,  an  additional 
simulation  including  X-rays  is  also  analyzed.  In  Figure  18,  the  inclusion  of  X-rays  results  in  a 
O.OlMeV  increase  in  average  photon  energy.  After  using  equation  (20),  total  activity  shows  an 
increase  to  605  Photon  Megacuries. 


Photon  Source  Definition  -  generated  from  ORIGEN  F71  file 
Filename:  today. f 71 

case  106.  3.600E03  seconds  0.000E00 

total  strength:  2.2393E19  gammas /second 
average  energy:  0.8220924456853023 

discrete  lines:  0.0000E00  gammas/second 
0.00%  of  energy 
average  line  energy:  0.0 

multigroup  bins:  2.2393E19  gammas/second 
100.00%  of  energy 

average  group  energy:  0.8220924456853023 

not  counted:  5. 5186E16  gammas/second 

1.10%  of  tne  original  energy 


Figure  18.  ORIGEN  Fallout  Analysis  Tool  output  for  average  energy  and  total 
activity  at  1  hour  for  a  90/10  U-235/U-238  fission  including  X-rays. 
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Because  the  SNC  is  independent  of  activity  density,  the  change  in  average  energy  and 
group  weighting  are  the  only  possible  impacts  of  X-ray  inclusion.  Nonetheless,  all  weighting 
values  used  in  the  calculation  of  SNCs  from  this  point  forward  combine  both  gamma  ray  and 
X-ray  contributions. 

SNC  for  90/10  U-235  Fission  Using  Multi-Group  Approach 

In  order  to  gauge  the  effect  of  the  energy  dependent  BUF  value  on  the  SNC,  the 
energy  spectmm  from  O.OMeV  to  3.5MeV  is  broken  up  into  35  energy  groups  (see  Figure  19). 
In  order  to  gain  increased  resolution  over  the  region  where  the  BUF  varies  significantly,  the 
range  of  O.OMeV  to  O.lOMeV  is  further  divided  into  10  subgroups. 


Figure  19.  ORIGEN  Fallout  Analysis  Tool  Photon  Source  output  showing  the 
distribution  of  energy  and  intensity. 
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Using  the  process  described  in  Chapters  II  and  III,  the  BUF  and  SNC  values  are 


calculated  for  each  individual  energy  group  (see  Figure  20)“. 


Once  calculated,  these  values  are  paired  with  weight  values  generated  by  the  ORIGEN  Fallout 
Analysis  Tool  and  combined  to  generate  a  total  SNC.  The  result  of  using  a  multi-group 
approach  to  calculate  a  SNC  for  the  spectmm  of  O.OMeV  -  3.5MeVat  a  time  of  one  hour  post 
detonation  is23: 


Multi-Group  SNC  =  8.6 x  10  16 


J-m2 
kg  -  sec-  Bq 


22  BUF  values  for  each  energy  group  are  shown  in  Figure  10 

23  Appendix  D  contains  all  tabled  weight  values  used  to  calculate  this  value 


36 


This  value  is  15%  greater  than  the  SNC  value  determined  in  Chapter  III.  The  source  of  this 


difference  is  discussed  in  “Addressing  the  Difference  in  Source  Normalization  Constant” 
section. 

SNC  for  90/10  U-235  Fission  Using  Average  Energy  and  BUF  Values 

After  determining  the  total  SNC  based  on  a  multi-group  approach  to  the  energy 
spectrum,  another  SNC  is  calculated  using  the  average  energy  of  0.81MeV  and  an  average  BUF 
value  for  the  energy  distribution.  The  average  BUF  value  is  calculated  using  the  weight  values 
generated  by  the  ORIGEN  Fallout  Analysis  Tool  combined  with  a  specific  BUF  value 
calculated  for  each  individual  energy  group.  This  yields  an  average  BUF  value  of  1.1. 24  The 
SNC  calculated  using  this  method  is25: 


Average  BUF  SNC  at  0.81MeV  =  8.6 xlCT16 


T  2 

J  -m 

kg-sec-Bq 


This  result  is  identical  to  that  of  the  multi-group  approach  in  the  previous  section.  This 
confirms  that  the  use  of  a  single  average  photon  energy  and  average  BUF  value  is  an  accurate 
approximation  in  calculating  the  SNC. 

Addressing  the  Difference  in  Source  Normalization  Constant 

Recall  the  SNC  calculated  in  Chapter  III  uses  an  average  gamma  ray  energy  of 
0.70MeV.  This  is  in  immediate  disagreement  with  the  ORIGEN  Fallout  Analysis  Tool 


24  Appendix  D  contains  all  tabled  weight  values  used  to  calculate  this  value 

25  This  calculation  shown  in  Appendix  A  (6) 
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estimate  of  0.81MeV.  If  the  method  of  calculating  the  SNC  contained  in  this  chapter  is 
accurate,  the  SNC  calculated  using  an  average  gamma  ray  energy  of  0.70MeV  should  match 
the  value  calculated  in  Chapter  III.  Therefore,  the  entire  process  from  the  previous  section  is 
repeated  using  an  average  energy  of  0.70MeV  and  average  BUF  of  1.1.  Using  these 
parameters,  the  SNC  is  determined  to  be26: 

f  J  —  m 2 

Average  BUF  SNC  at  0.70MeV  =  7.5  x  10  16  - 

kg  -  sec-  Bq 

This  value  represents  a  match  with  the  SNC  value  obtained  in  Chapter  III  and 
demonstrates  that  the  difference  in  SNC  values  are  due  to  the  use  of  0.81MeV  as  the  average 
photon  energy. 

Fallout  Distribution  Analysis  SNC  Results 

From  these  results,  two  observations  are  made:  First,  the  average  energy  of  a  fallout 
field  is  0.81MeV  instead  of  0.70MeV  as  calculated  by  the  ORIGEN  Fallout  Analysis  Tool. 
Second,  the  assumption  of  a  single  average  gamma  ray  energy  with  an  average  BUF  is  shown 
to  be  a  reasonable  approximation.  Although  the  BUF  changes  significantly  in  the  energy  range 
O.OMeV  to  0.50MeV,  the  relative  contributions  of  those  energy  groups  are  minor  in 
comparison  to  the  spectmm  as  a  whole.  However,  these  results  are  still  limited  by  the 
assumptions  made  in  the  development  of  this  model  (discussed  in  Chapters  II  and  III).  With 

26  Since  the  BUF  value  of  1.1  is  confirmed  in  this  chapter,  this  calculation  contains  identical  values  to  those  used  in  Chapter  III 
(page  21) 
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these  modifications  originating  from  the  fallout  distribution  analysis  in  this  chapter,  the  SNC  is 
further  refined  to  a  value  of: 


fission  yield  changes  depending  on  the  fuel  used.  This  fact  strongly  reinforces  the  use  of mks  units  as  the 
standard  in  discussing  SNCs. 
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V.  Varying  Activity  Density 

In  the  previous  chapters,  assumptions  regarding  the  source  normalization  constant 


(SNC)  were  analyzed.  As  shown  in  equation  (6),  activity  density  is  also  an  influential  term  in 
determining  dose  rate.  Stated  in  assumption  #4,  the  first  order  model  (FOM)  assumes  the 
activity  density  directly  beneath  the  detector  is  uniform  out  to  infinity.  Removing  the  spatial 
dependence  of  activity  density  allowed  it  to  be  removed  from  the  integral  as  a  constant: 


D[X,Y,l]  =  A[X,Y,l] 


Mjhvavg] 

V.  p  ) 


100  — 

r  e 


-ds 


(3) 


In  reality,  there  is  a  spatial  dependence  on  activity  density;  therefore,  this  assumption  must  be 
analyzed  to  gauge  its  impact  on  calculating  dose  rate. 

Activity  Deposition  on  the  Ground 

In  the  Introduction  to  the  Physics  of  Nuclear  Weapons  Effects,  Bridgman  discusses  the 
formulation  of  a  model  for  activity  distribution  on  the  ground  incorporating  both  vertical 
settling  and  horizontal  wind  transport  (Bridgman,  2001:41 1).  In  this  equation  the  lateral  spatial 
distribution  (modified  to  use  variables  consistent  with  this  thesis)  is27: 


A[r,t\  =  A[t] 


1 


y 
2l  a. 


42ft ( 


ft<J„ 


(20) 


Where, 


A[r,t]  represents  the  activity  density  at  a  position  r  [m]  and  at  time  t  [sec]  ; 

A[t ]  is  the  activity  density  at  time  t,  and 

cr.  is  the  standard  deviation  of  the  fallout  distribution  in  [m]. 


27  Bridgman  equation  (12-24) 
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This  expression  makes  several  assumptions.  First,  it  assumes  a  normalized  Gaussian 
distribution  in  horizontal  directions.  Although  toroidal  motion  during  cloud  rise  creates  a  ring 
shaped  distribution  (Jodoin,  1994:91-93),  ground  fallout  patterns  suggest  a  transition  to  a 
Gaussian  distribution  occurs  sometime  during  cloud  fall  (Bridgman,  2001:407).  For  purposes 
of  this  discussion,  we  also  assume  no  wind.  The  result  of  this  assumption  would  be  a  circular 
fallout  pattern  centered  on  ground  zero.  Although  this  assumption  is  unrealistic,  it  actually 
represents  the  worst  case  scenario  in  terms  of  activity  density  varying  with  respect  to  distance 
from  the  detector.  In  other  words,  a  “smeared”  Gaussian  distribution  only  in  the  horizontal 
axes  presents  a  larger  standard  deviation.  Finally,  we  assume  the  standard  deviation  of  fallout 
is  symmetric  —  further  leveraging  the  no  wind  assumption. 

The  distribution  represented  in  equation  (20)  is  inserted  into  the  integral  in  equation 

(3).  The  result  is  equation  (21),  where  the  normalization  factor  , — —  is  removed  to  create  a 

-> j27T(Jr 

correction  factor  that  maintains  a  value  between  0  and  1.  This  result  is  shown  in  Figure  21. 


Figure  21.  The  Gaussian  Activity  Correction  Factor 
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Thus  we  have  defined  a  correction  factor: 


1  5  -A 


rS  2 


ds 


Correction  Factor 


ds 


When  the  correction  factor  is  incorporated  with  equation  (3): 

-(ft[  hvm 


1  V-/r 


D[r,  1]  =  A[l\ 


f  P<XhvaVg^ 


V  P  Jai, 


1WJ 

r  e 


-t/.v 


(21) 


The  next  step  is  calculating  a  reasonable  value  for  cr . .  The  following  empirical 


relationship  (modified  to  use  variables  consistent  with  this  thesis)  relates  yield  to  cr. 


(Pugh, 1959:24): 


cr.  =  1.609c1 


0.70+1  in(  Y) 


3.25 

4.0+[ln(r)+5.4]2 


Where, 

cr r  represents  the  standard  deviation  of  the  fallout  distribution  in  [km],  and 
Y  is  yield  in  [kilotons]. 

Using  the  model  values  of  one  kiloton  and  a  detector  height  of  one  meter: 

cr.  =  2.9  km  =  2900 m 


(22) 


Non-Uniform  Activity  Results 

Evaluating  equation  (22)  using  an  average  photon  energy  of  0.81MeV  and  including  the  BUF 
value  of  1.1  obtained  in  Chapter  III,  the  SNC  is  calculated  to  be: 
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Variable  Activity  SNC  =  8.6x10 


-16 


J-m2 


kg-sec-Bq 


Calculated  in  Appendix  A.,  this  value  indicates  a  mere  1 .8x10  ~  %  decline  in  the  SNC  calculated 
in  Chapter  IV  from  varying  activity28. 

The  standard  deviation  of  the  fallout  distribution  is  the  driving  factor  in  any  impact  on 
the  SNC.  In  other  words,  the  smaller  the  standard  spatial  deviation  of  the  fallout  distribution, 
the  steeper  the  decline  in  activity  density  as  position  is  increased  from  the  detector.  The 
assumption  of  no  wind  represented  the  worst  case  scenario  (or  smallest  standard  deviation). 
This  means  at  worst  the  variation  in  activity  density  has  insignificant  impact.  Nonetheless, 
several  calculations  were  made  using  a  variety  of  values  for  standard  deviation  to  further 
illustrate  the  impact  on  correction  factor  Figure  (21)  and  SNC  (Figure  22). 


Figure  22.  Source  Normalization  Constant  as  a  function  of  fallout  standard 
deviation. 


28  See  Appendix  A  (7) 
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Shown  in  Figures  21  and  22,  even  as  standard  deviation  is  reduced  below  1km  the  correction 
factor  still  maintains  a  value  greater  than  0.99.  Next,  the  impacts  of  time  on  SNC  are 
considered. 
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VI.  The  Impact  of  Time  on  the  Source  Normalization  Constant 


As  shown  in  Chapter  III,  the  BUF  has  a  strong  energy  dependence.  In  Chapter  IV, 
the  reduction  of  the  fallout  energy  distribution  to  a  single  average  photon  energy  and  an 
average  BUF  is  shown  to  be  accurate  in  calculating  the  source  normalization  constant  (SNC). 
Though  as  time  progresses,  the  distribution  of  energy  is  expected  to  change  as  short-lived  high 
energy  emitting  isotopes  decay  into  longer-lived  lower  energy  emitting  isotopes.  This  change 
in  distribution  will  also  affect  the  average  BUF  value.  In  this  chapter,  we  investigate  the 
impacts  of  time  on  calculating  the  source  normalization  constant. 

The  Way-Wigner  Approximation 

In  1948,  Katherine  Way  and  Eugene  Wigner  published  what  is  known  as  the  Way- 
Wigner  approximation.  In  their  publication,  an  extensive  theoretical  treatment  of  radioactive 

decay  yields  the  relationship  (Way,  1948): 

DGD(t)  =  DGD(\yL2  (23) 


Where, 


DGD(t )  is  the  dose  rate  at  time  /in 


J 


kg  -  sec 


DGn  ( 1 )  is  the  dose  rate  at  unit  time  (lhr,  1  day,  1  month)  in 


J 


kg  -  sec 


and 


t  is  time  in  whatever  units  are  being  used  for  the  unit  time  value. 

Combining  equation  (23)  with  an  updated  version  of  equation  (6)  featuring  the  new  SNC  at 
one  hour  yields  equation  (24): 
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D[X,  Y,t\  =  8.6x1 0  16  A[X,  Y,  1]  ru  { - 

[  kg  -  sec 

This  relationship  reveals  an  important  distinction  in  the  application  of  the  Way-Wigner 
approximation  which  is  discussed  in  the  next  section. 

Time  Impact  on  Source  Normalization  Constant  (1  -  72  Hours) 

As  shown  in  equation  (24),  dose  rate  is  calculated  by  multiplying  an  activity  density  by 
a  source  normalization  constant.  The  cumulative  impact  of  time  affects  both  of  these  factors. 
Since  half-life  generally  increases  as  fallout  progresses  down  the  decay  chain,  activity  will 
inherently  decline.  In  fact,  this  is  precisely  the  effect  addressed  by  the  Way-Wigner 
approximation.  Unfortunately,  this  does  not  address  any  impacts  on  the  source  normalization 
constant. 

For  example,  as  the  activity  density  changes  the  distribution  of  that  activity  across  the 

energy  spectmm  also  changes.  This  means  the  average  photon  energy  for  the  spectrum  will 

change.  As  a  consequence,  this  also  changes  the  weight  values  used  to  calculate  an  average 

BUF  for  the  distribution.  The  question  becomes:  are  these  changes  in  distribution  and  energy 

significant  enough  to  warrant  a  separate  time  correction  term? 

To  determine  the  impact  of  time  on  the  SNC,  runs  for  90/10  U-235/U-238  at  various 
times  are  created  within  the  ORIGEN  Fallout  Analysis  Tool  for  times  of  12  hours,  24  hours, 
48  hours,  and  72  hours.  Identical  to  the  procedure  used  in  Chapter  IV,  each  time  interval  is 
analyzed  to  obtain  average  photon  energy  and  a  weighted  distribution  function  to  calculate  an 
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average  BUF.  A  summary  of  these  results  is  illustrated  in  Figure  24.  A  sample  of  completed 
code  results  is  included  in  Appendix  D  for  the  12  hour  case. 


Figure  24.  The  source  normalization  constant  as  a  function  of  time 


As  shown  in  Figure  24,  the  SNC  decreases  steadily  over  time.  As  time  increases,  the 
declining  average  photon  energy  and  the  impact  of  different  weight  distributions  on  the 
average  BUF  are  the  most  likely  causes;  however,  further  analysis  is  necessary  to  determine 
which  of  these  factors  drives  the  decline. 

If  fallout  activity  distributions  increasingly  favor  low  energies  over  time,  the  average 
BUF  value  is  expected  to  increase  (see  Figure  10).  This  prediction  holds  true  when  average 
BUF  values  are  calculated  with  weighting  factors  from  the  ORIGEN  Fallout  Analysis  Tool 
(Figure  25).  An  important  observation  is  the  range  of  increase  beginning  at  1.1  and  rising  to  a 
value  of  1.2  before  72  hours,  representing  a  9.0%  increase. 
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Average  BUF  vs.  Time 
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Figure  25.  The  build  up  factor  as  a  function  of  time 


As  a  consequence  of  the  energy  distribution  initially  favoring  lower  energies  as  time 
increases,  the  average  photon  energy  is  also  expected  to  decline.  This  is  verified  again  by 
graphing  results  obtained  from  the  ORIGEN  Fallout  Analysis  Tool  (Figure  26). 


Figure  26.  Average  photon  energy  as  a  function  of  time 
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In  examining  Figure  26,  a  decline  in  average  energy  from  0.81MeV  to  0.45MeV  is 
observed  over  the  first  72  hours.  This  represents  a  56%  decline,  overcoming  the  effects  of 
increasing  BUF  values  discussed  earlier. 

Time  Impact  on  Source  Normalization  Constant  (72  Hours  -  1  Month) 

Figures  24-26  seem  to  indicate  an  approach  to  stable  values  after  72  hours;  however, 
additional  analysis  of  average  BUF  and  average  energy  extended  out  to  one  month  indicates 
differendy.  Using  the  same  methods  as  previously  described  for  90/10  U-235/U-238  fission, 
Figures  27-29  show  a  steady  increase  in  average  photon  energy,  a  decrease  in  average  BUF, 
and  an  increase  of  the  SNC  value: 


Average  Photon  Energy  (MeV) 

ppppppppo 

DPMUI^^OlsJOoiD 

_ ▲ _ 

Average  Photon  Energy  vs.  Time 

0  100  200  300  400  500  600  700  800 

Time  (hours) 

Figure  27.  Average  Photon  Energy  as  a  function  of  time  out  to  one  month. 
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Average  BUF  vs.  Time 
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Figure  28.  Average  Photon  Energy  as  a  function  of  time  out  to  one  month. 


Using  previously  described  methods,  the  SNC  at  one  month  is: 


SNC  at  1  Month  =  7.4x10  16 


f  ,  2  1 

J  -m 

< - > 

kg  -  sec-  Bq 
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Figures  27-29  reinforce  the  fact  that  average  photon  energy  is  the  driving  factor  in  determining 
the  SNC,  but  what  causes  the  average  energy  to  increase  after  72  hours? 

Although  the  entire  energy  distribution  varies  over  time,  one  isotope  in  particular  is 
determined  to  be  a  major  contributor  to  the  increase  in  both  average  gamma  ray  energy  and 
SNC.  Consider  the  decay  chain: 

140 Xe  -»  140 Cs  ->  140 Ba  ->  140 La 

Half-lives  for  the  first  three  isotopes  are  13.6  sec,  1.06  min,  12.75  days,  respectively.  The  half- 
life  of  140  La  is  only  1.68  days,  meaning  140  La  and  140  Ba  will  form  a  transient  equilibrium. 
Furthermore,  140  La  emits  a  1.60MeV  gamma  ray  compared  to  a  0.54MeV  gamma  ray  from 
140 Ba  (Baum,  2002:56). 

According  to  the  ORIGEN  Fallout  Analysis  Tool,  the  energy  contribution  of  14uLa  at 
a  time  of  72  hours  is  2.6%  of  the  entire  fallout  energy  distribution^.  As  time  increases,  the 

contribution  of  140  La  also  continues  to  increase.  At  two  weeks,  14uLa  accounts  for  13%.  At 
one  month,  140  La  accounts  for  16%.  Although  several  other  decay  chains  influence  the 
average  gamma  ray  energy,  these  results  suggest  140La  is  a  major  contributor.  Using  Figure  27 
and  29  as  a  reference,  the  initial  decline  in  average  gamma  ray  energy  (and  SNC)  can  be  viewed 
as  a  lack  of  140  La  present.  After  a  few  days,  140  Ba  begins  to  decay  into 140  La  .  The 
comparatively  short  half-life  of  1.68  days  coupled  with  the  1.6MeV  gamma  ray  emission  of 
140  La  pushes  the  average  gamma  ray  energy  back  up  to  0.65MeV,  thereby  increasing  the  SNC 
to  the  aforementioned  value. 

29  See  Appendix  C 
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VII.  Conclusions  and  Recommendations 


As  previously  stated,  the  objective  of  this  research  was  to  analyze  the  assumptions 
listed  in  Chapter  I  regarding  the  first  order  model  developed  to  calculate  dose  rate  of  a 
detector  positioned  one  meter  above  a  fallout  field.  As  the  assumptions  were  relaxed,  any 
quantitative  impacts  were  combined  to  form  an  updated  source  normalization  constant. 

In  Chapters  II  and  III,  assumption  #1  of  “no  scattering”  was  partially  relaxed  by 
utilizing  an  integral  method  of  successive  scatters.  By  including  one  scatter  and  a  limited 
second  scatter  approximation,  the  build  up  factor  (BUF)  was  determined  to  be  1.1  for  an 
average  photon  energy  of  0.70MeV.  This  result  indicates  scattering  contributes  an  additional 
10%  to  dose  rate. 

In  analyzing  BUF  values  for  energies  other  than  0.70MeV,  it  became  apparent  that  the 
BUF  changes  significantly  over  the  energy  spectmm  -  most  notably  to  a  value  of  1.9  at 
0.02MeV.  Since  the  location  of  this  BUF  peak  is  an  artifact  of  only  considering  a  single 
scatter,  further  research  into  multiple  scatters  is  highly  recommended.  A  Monte-Carlo  scheme 
would  be  best  suited  in  this  endeavor. 

Since  the  BUF  was  shown  to  carry  significant  energy  dependence,  the  assumption  of 
using  a  single  average  photon  energy  to  represent  the  full  fallout  energy  distribution  was 
addressed  next.  In  Chapter  IV,  individual  BUF  values  were  assigned  to  energy  groups 
spanning  the  range  of  O.OMeV  to  3.5MeV.  Using  the  ORIGEN  Fallout  Analysis  Tool 
provided  by  Oak  Ridge  National  Laboratory,  weighted  values  for  the  fallout  energy 
distribution  were  obtained  for  each  energy  group  and  used  to  calculate  a  weighted  BUF 
average.  The  result  of  this  effort  was  a  confirmation  of  a  BUF  value  of  1.1. 
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Over  the  course  of  determining  an  average  BUF  value,  it  was  shown  by  the  ORIGEN 
Fallout  Analysis  Tool  that  the  average  photon  energy  of  a  fallout  field  at  one  hour  is  0.81MeV 
-  in  contrast  with  the  widely  accepted  value  of  0.70MeV  listed  in  The  Effects  of  Nuclear  Weapons 
(Glasstone,  1977:454).  Using  the  new  average  photon  energy  and  average  BUF  value,  a 
refined  source  normalization  constant  at  one  hour  was  determined  to  be: 


In  comparing  this  SNC  with  one  calculated  using  a  multi-group  approach,  it  was 
determined  that  fallout  distributions  can  be  modeled  using  a  single  average  photon  energy  in 
conjunction  with  a  weighted  average  BUF  term.  This  was  proven  by  identical  SNC  values 
calculated  with  the  two  methods. 

Chapter  V  explored  the  impact  of  non-uniform  activity  across  a  fallout  field.  Using  a 
normalized  Gaussian  distribution  to  model  a  no  wind  (worst  case)  scenario,  less  than  a 
1.8x10  2  %  decline  was  observed  on  the  SNC.  Further  analysis  of  standard  deviations  from 
1000m  to  10,000m  showed  the  impacts  of  standard  deviation  on  the  SNC  and  correction 
factor,  confirming  the  result  that  assuming  uniform  activity  density  (assumption  #4) 
determined  by  the  activity  density  directly  beneath  the  detector  yields  an  accurate 
approximation  of  SNC  within  a  fraction  of  a  percent. 

In  Chapter  VI,  the  effects  of  time  (assumption  #3)  were  explored  to  determine  if  a 
separate  time  correction  factor  was  needed  to  address  the  changes  in  energy  distributions  that 
compose  SNC  and  BUF  values.  After  analyzing  the  effects  of  time  on  SNC,  a  separate  time 
correction  factor  was  found  to  be  necessary  to  adjust  the  SNC  at  times  other  than  one  hour. 
This  was  determined  by  discovering  a  sudden  57%  decline  in  SNC  over  the  first  72  hours, 
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followed  by  an  increase  to  86%  of  its  original  value  after  one  month.  This  result  indicates 
assumption  #3  is  not  sufficient  in  modeling  the  effects  of  time  on  dose  rate  and  an  additional 
time  correction  factor  for  the  SNC  is  required. 

Furthermore,  140  La  is  determined  to  be  a  major  contributor  to  the  increase  in  average 
gamma  ray  energy  and  SNC  at  times  greater  than  72  hours.  Using  the  ORIGEN  Fallout 
Analysis  Tool,  the  energy  contribution  of  14<)La  is  shown  to  increase  from  2.8%  to  16%  of  the 
fallout  energy  distribution  by  the  end  of  the  first  month.  This  large  contribution  coupled  with 

the  1.6MeV  gamma  ray  emission  of  140  La  drives  the  average  gamma  ray  energy  and  SNC  to 
increase  after  72  hours. 

Finally,  the  assumption  of  ground  as  a  smooth  and  level  plane  has  already  been 
addressed  in  The  Self-Shielding  of  Fallout  Gamma  Rays  by  Terrain  Roughness  by  Herte.  Using  a 
Monte-Carlo  code  (MORSE)  to  model  surface  roughness,  his  results  indicate  self-shielding 
between  2-5%  for  soil.  This  value  is  shown  to  vary  significantly  based  on  surface  type,  so  its 
effects  are  not  incorporated  into  the  final  SNC  value. 

Although  several  aspects  regarding  the  SNC  were  analyzed,  several  assumptions  are 
left  outstanding.  As  previously  mentioned,  this  analysis  addresses  only  one  scatter  and  a 
limited  approximation  of  second  scatter.  Table  1  demonstrates  the  impact  of  this  limitation. 
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Table  1.  Flight-by-Flight  Average  Photon  Energies 


Flight 

Incident  Energy  (MeV) 

Average  Post  Scatter  energy  (MeV) 

Mean  Free  Path  in  Air  (m) 

0 

0.70 

0.38 

109 

■ 

0.38 

0.25 

83 

0.25 

0.18 

72 

0.18 

0.14 

67 

0.14 

0.11 

59 

0.11 

0.093 

57 

0.093 

0.078 

52 

7 

0.078 

0.070 

48 

8 

0.070 

0.062 
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As  shown  in  Table  1,  average  photon  energy  decreases  an  order  of  magnitude  after  8 
scatters,  but  mean  free  path  is  only  reduced  42%.  This  result  strongly  indicates  that  modeling 
many  scatters  is  necessary  to  refine  the  BUF  term.  Furthermore,  most  of  the  analyzed  scatters 
were  assumed  to  be  isotropic  -  clearly  disproven  by  Klein-Nishina.  While  certain 
approximations  made  rough  adjustments  to  correct  for  non-isotropic  scatter,  a  more  in-depth 
analysis  may  show  further  impact  to  the  BUF  term. 


55 


Appendix  A 
Supporting  calculations 


(1)  Post-Scatter  energy  for  0.70MeV: 

Using  equations  from  The  Atomic  Nucleus  (Evans,  675),  the  relationship  between 

incident  photon  energy,  scattering  angle,  and  post-scatter  energy  is  defined: 

hv  1 

hv0  l  +  a(l-cos.9) 


As  discussed  in  Chapter  II  a  weighted  average  is  used  to  determine  the  average  scattering  angle 

3. 
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dS 

=  1.17  radians 


This  average  scattering  angle  is  then  used  in  equation  (16)  to  find  average  post-scatter  energy: 


hv'  = - ^hvn  =  0.38  MeV 

l  +  a(l-cos5) 


The  Mathematica  script  for  this  calculation  is: 

(2.81794*  10~15)2 


r  i 

NIntegrate  - 


(ER[S])J«  (ER [(5] +ER[r>]  -Sin[(J]J)  2  7r.Sin[5]  *<J,  {<?,  0,  7t}]  / 


r  (2. 81794  .lO'15) 2  ,  ,  ,  , 

NIntegrate  |  - - —  »  (ER[l?])2  *  (ER[tf]'1  *  ER[0]  -  Sin  [Up)  2  7r*Sin[<J]  ,  {6,  0,  7t}J 


1.17366 


1 

- *hv0 

1  +  a*  (1  -  Cos  [%] ) 

0.380429 
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(2)  Second  Flight  (1st  scatter)  contribution  to  the  detector  of  7.5%: 

Several  constants  are  required  for  this  calculation,  starting  with  the  incident  energy  and  average 
first/ second  scatter  energies: 

hv0  =  O.lOMeV 
hv'  =  0.38  MeV 
hv”  =  0.25  MeV 


The  cross  sections  used  for  this  calculation  are: 

'  n°  =9.2xl0~3{/tr1} 

i/'  =1.2x1 0  2 1 m  1 1 

A  =13{m-1} 

/u's  =\1  {rn1} 

H°a  =  7.1x10  3  {m 

ju'a=  3.5x1  O’2  {m1} 

/u°  =13{w-1} 

//;  =  17  {m-1} 

//;=1.2xiU‘|m  ) 

//"=  1.4x1 0"2  {m-1} 

As  previously  stated,  density  for  ground  and  air  are  taken  as: 

Air:  p  =  1.23  j-^yj 
Ground:  p  =  1900  j^-j 

Evaluating  equation  (12)  from  Chapter  II  yields  a  value  of  isotropic  air  scatter  contribution  to 
the  detector: 


For  Air:  < 


'  *  t 

H"s  =  1.4x1 0^2{m" 
p°a  =  9.9x  10“7  jm-1 
p'a  =  4.9x1  O’6  {m'1 
//''  =  1.7xl0'5{m  1 
p°  =9.2x1 0“3{m_1 

,  .  »  _ _ _i  i 


For  Ground:  i 
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DZ=C 


n  'jo 

J  Ex  ju°  (h  -«)]  ZT,  [  jLi'u  ]  du  +  J  [  H°  ( w  +  h)]px  [  //'vv]  dw 


=  4,5.9x10  17 


J 


kg- 


1  /  (AMuAl  +  cs[hvl])  \ 

DoselSA  =  -  *  hvl  «  - : -  *  Agd  .  MeVtoJ  *  AMuSO  . 

4  \  ARoe  J 

(HIntegrate  [F [AMuTO  (1  -  x)  ]  *F[AMuTl*  (x)  ]  ,  {x,  0,  1},  WorkingPrecision  -.  10]  + 

NIntegrate  [F[  AMuTO  (w  +  1)  ]  *  F [AMuTl  *  (w)  ] ,  {w,  0,  Infinity},  WorkingPrecision-*  10]) 

5.85205  x  10'17 

Likewise  for  ground  contribution,  evaluating  equation  (1 5)  yields  a  value  of  isotropic  ground 
scatter  contribution  to  the  detector: 


Dgd  =i 

^ SF 


cSE<  Mgd  ("-*> 


Mi 


gd 


(w -/z)  +  (///  )  h 

\  /  air 


du  =  Agd  4.8x10  17 


J 


kg  -  sec 


1  /  (AMuAl  +  ca[hvl])  \ 

DoselSG  =  -  *  MeVtoJ  «  hvl  * - -  »  Agd  »  GMuSO  * 

4  1  ARoe  I 

HIntegrate  [F[GMuT0  (x-1)]  *  F[  (GHuTl  *  (x  -  1) )  +  (AMuTl)],  {x,  1,  Infinity),  WorkingPrecision-.  10]  | 
4 . 80302  x  10"17 


Combined  with  the  evaluated  value  of  direct  flight  (DF)  contribution  in  equation  (4): 


DDF=Agd  6.8xl016 


J 


kg  -  sec 


/ 1 

DoseDF  =  -  *  hvO  »  MeVtoJ  «  Agd  *  - 
2  \ 

6.7S885xlO'lc 


ARoe 


«  F [AMuTO] 


The  total  dose  rate  is  the  sum  of  these  three  values.  Various  percentage  based  contributions 
can  be  calculated  from  these  values.  For  example: 
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SF 


b£+bfF+bDF 


=  %  Air  Scatter  Contribution  =  7.5% 
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=  %  Ground  Scatter  Contribution  =  6.1% 


sec 
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DoselSA 


.100 


DoselSG 


.100 


(DoselSA  +  DoselSG  +  DoseDF)  (DoselSA  +  DoselSG  +  DoseDF) 

7.47927  6.13855 


(3)  The  contribution  of  air  scatter  coming  from  the  region  above  the  detector: 

Evaluating  equations  (10)  and  (1 1)  yield  the  dose  rate  contributions  from  air  scatter  received 
from  below  and  above  the  detector  (respectively).  When  combined,  these  regions  become  the 
total  contribution  of  air  scatter  to  the  detector  equation  (12).  Using  the  constant  value 
obtained  from  equation  (12)  as  the  total  dose  rate,  the  percentage  contributions  from  the 
region  above  or  below  the  detector  can  be  calculated. 
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00  ')airs"  3 
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For  Region  III  (Above):  Z)“"  (///)  =  C 


}  } - 7 - ds' +  J  - r — ds" 
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h'—h 


dti 


Combine  Region  II  and  III  for  total  air  contribution: 


Dmr  =  C 


n 

J  iq  [ pi’ ( h  - u)~^E1  [ pi'u]  du  +  J  iq  [ Et  (w  +  h)~^El  \ /C vv]  dw 


The  percent  contribution  coming  from  Region  III  (Above)  is  then: 


Dmr 

^ SF(IU ) 

7~)air 

USF 


=  83% 


NIntegrate  [F  [AMuTO  (w  +  1)  ]  .  F  [AMuTl  (v)  ]  ,  {w,  0 ,  Infinity}  ,  WorkingPrecision  -*  10]  / 

( NIntegrate [ F [ AMuTO  (1-x)]  .F[AMuTl  (x)],  (x,  0,  1},  WorkingPrecision  10]  + 

NIntegrate  [F  [AMuTO  (w  +  1)  ]  *F[AMuTl  (w)  ]  ,  {w,  0,  Infinity},  WorkingPrecision-.  10])  .100 

82.9597249 
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(4)  The  BUF  of  1.1  for  average  incident  energy  of  0.70MeV: 

This  value  was  obtained  using  equation  (1 8)  of  Chapter  2.  As  shown  below,  equation  (1 8)  was 
modified  to  include  the  Ground  Non-isotropic  Adjustment  Term. 


3  K 


2  2 
fZL 

'  hv' ' 

3 

'hv,' 

+ 

"A 

~s; 

i  2 

KhV0j 

_yhV  ) 

[hvj 

GNAT  = 


2 


In  2 

f  !o_ 

'  hv'  > 

3 

s~ 

J  2 

0 

yhvoJ 

_  V  hv'  J 

f  hA 
\  K  j 


d& 


d& 


o.io 


r  (2.81794  „10-15) 2  3,1  2,  (  n  3*11  / 

GHAT  =  NIntegrate |  — - - - —  *  (ER[  9])  »  (ERfO]'1  +  ER['3]  -  Sin[(S]2) ,  |<5,  -  ,  —  JJ 

r  (2. 81794. 10-15) 2  ,  ,  ,  , 

NIntegrate  y  — - - - —  «  (ER  [<?] )  *  (ER  [$] _1  +  ER  [t?]  -Sin[(J]‘),  [&,  0,  2  7r}  J 

0.0992015 


This  term  was  used  to  determine  the  fraction  for  ground  scatter  near  the  detector  that  was  able 
to  scatter  back  into  the  air.  From  this  value,  a  limited  second  scatter  approximation  (LSSA) 
was  made. 


DLSSA=GNATxhv" 

(K\ 

4 

l  P) 

Assa=^A  9x10  18  ■ 

J 

4 A4)a 


kg  -  sec 


1 


00 


J  E\  [4(1  -  «)]  E\  [,ii”(u  )]du  +  J  A,  [//,'( w  + 1)]  A,  [it"(w)]dw 

0  0 


1  /  (AMuA2  +  [hv2] )  \ 

DoselSAlSG  =  GNAT  *  -  *  MeVtoJ  *  hv2  *  -  »  Agd  *  AMuSl  * 

4  \  ARoe  I 

(NIntegrate [F[AMuTl  (1  -  x)  ]  «F(AMuT2  (x)  ]  ,  {x,  0,  1},  WorkingPrecision 6]  + 

NIntegrate  [F[AMuTl  (w  +  1)  ]  *F[AMuT2  (w)  ] ,  (w,  0,  Infinity},  WorxingPrecision  -*  6]) 


3 . 87647  x  10”ls 


Combining  these  values  with  previously  calculated  values  for  air  and  ground  contributions,  the 
BUF  is  calculated  using  equation  (18): 


60 


BUF  = 


1 1  j;;+(cjv/irxDg)+4m  t  ] 


(5)  Calculating  Mean  Free  Path: 

Calculating  the  mean  free  path  for  a  photon  is  simply  the  inverse  of  its  macroscopic  cross 
section  for  attenuation.  For  a  0.02MeV  gamma  ray  in  air: 


Ht  =5.5x10  2 1  m 


—  =  1 8  \m 
F, 


For  a  0.70MeV  gamma  ray  in  air: 


//,  =  9.2x1 0“3  {/»-*} 


—  =  \09\m 

Ft 


For  a  0.02MeV  gamma  ray  in  ground: 


//,  =  475.5{m _1} 
1 


=  2.1x10  3U 


Ft 


For  a  0.70MeV  gamma  ray  in  ground: 


jut  =12.8{m_1} 

—  =  7.8x10  2 
F, 


(6)  Calculations  involving  average  energy  of  0.81MeV: 

Several  constants  are  required  for  this  calculation,  starting  with  the  incident  energy  and 
average  first/ second  scatter  energies: 
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hv0  =  OMMeV 
hv'  =  0A\3MeV 
hv"  =  0.261  MeV 


The  cross  sections  used  for  this  calculation  are: 


For  Air:  < 


//;  =8.661x10  3  {/it 
ju's=  1.1 54x1 0"2  [m "'} 
//;  =  1. 365x1  O'2  {777-1} 
/u°  =  7.167x10  7  {/77  '} 
jii'a  =  3.850 x  10  6  {777  '} 
//;  =  1.415x10  5  {777^} 
n°  =  8.64x  10  3  {777  '} 
//;  =  i.i5xio 2 {777 '} 
//,"=  1.36x10"  {777  '} 


For  Ground: 


//;  =  12.02(777-'} 
ju's  =  16.01(777-'} 
p°a  =  5.1 16xl0“3  {777  '} 
p'a  =2.752xl0-2{  777-'} 

n°  =12.02(777-'} 
n't=  16.05(777-'} 


As  previously  stated,  density  for  ground  and  air  are  taken  as: 

Air:^’-230) 

Ground:  p  =  190o|— ^-| 


Evaluating  equation  (12)  from  Chapter  II  yields  a  value  of  isotropic  air  scatter  contribution  to 
the  detector: 


DZ=C 


n 

| Ex  (h  -  «)]  El  [ pit ]  du  +  J El  [} ju° (vv  +  /7)>,  [ p'tw\  dw 


=  Agd  6.2xl017- 


kg 


1  /  (AMuAl  +  oa  [hvO] )  \ 

DoselSA  =  -  .  hvl »  -  *  Agd  *  MeVtoJ  *  AMuSO  * 

4  1  ARoe  7 

(NIntegrate [F [AMuTO  (1  -  x)  ]  »F[AMuTl  (x)  ] ,  {x,  0,  1},  WortingPrecision 10]  + 

NIntegrate[F[A14uT0  (w  +  1)  ]  *  F  [  AHuTl  (w)  ] ,  {w,  0,  Infinity},  WortingPrecision-*  10]) 


6 . 15055  x  10-17 


sec 
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Likewise  for  ground  contribution,  evaluating  equation  (15)  yields  a  value  of  isotropic  ground 
scatter  contribution  to  the  detector: 


Dfr  =cf£, [(//,")  (u-K)\e\(A  (u-h)  +  (A 

i  L  £a  J  _  V  /  gd  V  /  air 


du  =  v4grf5.1xl0  17 


J 


kg  -  sec 


DoselSG  =  -  *  MeVtoJ  *  hvl  * 

4 


(J 


:) 


«  Agd  »  GMuSO  * 


ARoe 

NIntegrate [F [GMuTO  (x -  1) ]  »F[((34uTl«  (x-1))  +  (AMuTl) ]  ,  {x,  1,  Infinity},  WorfcingPrecision -►  10] 
5.07239  x  lO'17 


Combined  with  the  evaluated  value  of  direct  flight  (DF)  contribution  in  equation  (4): 

J 


DDF=Agdl.lxW 


-16 


kg  -  sec 


1  /  (AMuAO  +  a.rhvO])  \ 

DoseDF  =  -  .  hvO  «  HeVtoJ  «  Agd  *  -  *  F  [  AMuTO  ] 

2  \  ARoe  / 


7.73314x10' 


This  value  was  obtained  using  equation  (1 8)  of  Chapter  2.  As  shown  below,  equation  (1 8)  was 
modified  to  include  the  Ground  Non-isotropic  Adjustment  Term. 

3  7T 


2  2 

C- 

fhv^ 

3 

o| 

^  | 

+ 

f— 1 

i  2 

\hv  J 

J1V0, 

-  sin2  9 


d9 


GNAT  =  ■ 


In  2 


rhy^ 

o  -  \hvoJ 


k  lit,  \  f  fov'  4 


hv( 

\hv  j 


\  K  j 


-sin'  9 


=  0.10 


d9 


GNAT  =  NIntegrate 


(2. 81794. 10'15) 2 


.  (ER [i?] ) 3  .  (ERtd]'1  +  ER [d]  -  Sin[3]2) ,  {<5,  -  ,  — }]/ 


,  (2. 81794. 10'15} 2  ,  ,  , 

NIntegrate  [- - —  *  (ER  [0] ) 3  .  (ER  [6] +ER[t.»]  -  Sin[0]2) ,  {«,  0,  2  7t}j 


0.0875075 
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This  term  was  used  to  determine  the  fraction  for  ground  scatter  near  the  detector  that  was  able 
to  scatter  back  into  the  air.  From  this  value,  a  limited  second  scatter  approximation  (LSSA) 
was  made. 


DLSSA=GNATxhv" 


DLssa  —  ^gj3.9xl0 


P  Jai 

J 


AAP's)a 


1  W 

J Ex  [ju't(\  -u)\ E,  [jii''{u)\lu  +  J El  [ju't{w  + 1)] El  [u"(w)]clw 


kg  -  sec 


1  /  (A14uA2  +  a,  [Irol] )  \ 

DoselSAlSG  =  GNAT  .  -  *  HeVtoJ  .  hv2  .  - - -  »  Agd  .  AlluSl . 

4  1  ARoe  I 

( NIntegr ate  [  F  [  AHuTl  (1  -  x)  ]  *F[AMuT2  (x)  ] ,  {x,  0,  1},  WorkingPrecision  -♦  6]  + 

NIntegrate  [F  [AMuTl  (w  +  1)  ]  „F[AMuT2  (w)],  {w,  0,  Infinity},  WorkingPrecision  ->  6]) 

3.82763  x  10~ls 


Combining  these  values  with  previously  calculated  values  for  air  and  ground  contributions,  the 
BUF  is  calculated  using  equation  (18): 


BUF  =  l  + 


D‘sr  +  (GNAT  xD^.)  +  bISS4 


D 


=  1.1 


DF 


Using  the  average  BUF  and  direct  flight  dose  rate  calculated  above  (with  an  assumed  activity 

density  of  1  {Bq/m2}  to  yield  SNC)  results  in: 

SNC  =  DoseDF  *1.1 
8 . 59999  x  10‘le 


(7)  The  Impact  of  Varying  Activity  Density  on  SNC 

In  Chapter  V,  a  correction  factor  for  varying  activity  density  is  defined  as: 

if 

•o  2(  ar2  j 

- - - ds 

i  s 

Correction  Factor  =  - - - - - 

- ds 

J  C 


64 


Incorporating  this  correction  factor  into  equation  (3)  yields: 


1  S  -h 


f  U  [ llV 

D[r,\\  =  A[  1]  a  avg 


Pa[hvave]]  _  rre-(M^s)e-2{^ 


V  P  ) 


-ds 


1  s~-lr 


Where  e  v  r  7  is  the  correction  factor  for  varying  activity  density.  As  shown  in  previous 
calculations,  the  SNC  without  the  correction  factor  is: 


1 1 

SNC  =  -  *  hvO  *  MeVtoJ  »  Agd  *  [  - 
8 . 59999  x  10'16  Agd 


ARoe 


•) 


*  F[AMuT0]  »  BUF 


When  the  correction  factor  is  incorporated  into  this  function,  the  result  is: 

[$] 


-(AHaT0*s)-.5 


/  i 

SNC  =  -  *  hvO  *  MeVtoJ  *  Agd  *  [  - 


8.5984  xlO'16  Agd 


ARoe 


•)-r 


•  ds  //  N 


The  change  in  SNC  is  calculated  by: 


(1  -  (8.5984*A-16/8. 59999*A-16  ) )  *100 
0.0184884 


ASNC  =  1.8x10' 2% 
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Appendix  B 

ORIGEN  Fallout  Analysis  Tool  Parameters 

We  begin  with  opening  the  Fallout  Analysis  Tool.  The  User  Interface  (UI)  should 
look  like  Figure  30  below: 


\4>\  Fallout  Analysis  Tools 
File  Utilities  Help 

1.  Edit  Model  Parameters 

f  Edit  | 

2.  Run  Model  Using  SCALE 

|  Run  | 

3.  View  SCALE  Files 

(3)  ARP/ORIGEN-S  input  fHe 
SCALE  screen  messages 
SCALE  error  messages 
Full  text  output  file 
Explore  all  results 

4.  Advanced  Analysis 

I  Start  I 


Figure  30.  The  ORIGEN  Fallout  Analysis  Tool  User  Interface 


Under  “File”  on  the  toolbar,  select  “New.”  This  will  open  another  screen  allowing  the  user  to 
edit  various  parameters  about  the  model.  In  Figure  31,  the  parameters  for  the  90%/ 10%  U- 
235/U-238  fission  are  used. 
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Figure  31.  ORIGEN  Fallout  Analysis  Tool  Parameters 


As  shown  in  Figure  31,  a  1.0KT  yield  with  a  total  mass  of  1.0kg  is  used  in  this 

simulation.  The  values  for  duration,  usable  energy  per  fission,  steps  during  detonation  and 

decay  steps  per  decade  are  all  default  values,  but  may  be  altered  if  desired. 

The  “decay  times”  block  allows  the  user  to  input  various  times  of  interest  to  focus 

analysis.  Common  denominations  such  as  1  hour,  1  day,  or  1  week  may  be  added  using  the 

“Easy  Add  Time”  dropdown  menu.  Times  can  also  be  manually  added  (in  seconds)  using  the 

“Add”  button  on  the  bottom  of  the  block.  Once  all  parameters  are  saved,  the  user  is  returned 

to  the  main  UI.  From  the  Main  UI,  the  user  now  selects  the  “Run”  button  in  block  2.  After 

10-20  seconds,  new  options  will  become  available  under  block  3,  “View  Scale  Files. 

Select  “Explore  all  results”  and  click  view  at  the  bottom  of  the  block.  This  opens  an 

analysis  UI  with  several  useful  features.  For  purposes  of  this  research,  the  user  should  select 
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the  case  time  in  which  they  are  interested  (3600  seconds),  select  “Export”  from  the  toolbar, 
and  then  “Photon  Source.”  (Figure  32) 


Figure  32.  The  File  Explorer 


Once  selected,  the  export  photon  source  option  with  display  the  final  UI  screen 
allowing  the  user  to  specify  various  parameters.  The  parameters  used  for  the  research 
contained  within  this  document  are  illustrated  in  Figure  32.  Once  listed,  the  user  selects  the 
“compute”  button  at  the  bottom  of  the  interface. 
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Figure  33.  The  Export  Photon  Source  UI 


Once  computed,  a  plot  and  a  text  file  will  be  generated  and  may  be  saved  if  desired. 
As  shown  in  Figure  33,  the  average  energy  for  90/10  fission  at  a  time  of  1  hour  is  0.82MeV. 
Also  useful  is  the  “not  counted”  line  which  gives  the  total  activity  and  percent  of  original 
energy  not  contained  within  the  user  specified  energy  range.  For  example,  the  range  of  0.0  - 
3.5MeV  contains  roughly  99%  of  the  entire  spectrum. 
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Appendix  C 

ORIGEN  Fallout  Analysis  tool  Output 


100%  Pu-239  at  1  Hour 


Photon  Source  Definition  -  generated  from  ORIGEN  F71  file 
Filename:  today. f71 

case  106.  3.600E03  seconds  0.000E00 

total  strength:  2.2051E19  gammas/second 
average  energy:  0.7865806796307987 

discrete  lines:  0.0000E00  gammas/second 
0.00%  of  energy 
average  line  energy:  0.0 

multigroup  bins:  2.2051E19  gammas/second 
100.00%  of  energy 

average  group  energy:  0.7865806796307987 

not  counted:  1.4098E11  gammas/second 

0.00%  of  the  original  energy 


Multigroup  bins  (in  MeV)  and  their  probabilities 


groups:  50 

fraction : 

5 . 0000E00 

4 . 6645E-09 

4 . 9000E00 

4 . 1081E-07 

4 . 8000E00 

6. 5455E-06 

4 . 7000E00 

9. 4904E-09 

4 . 6000E00 

2 . 57  97E-08 

4 . 5000E00 

6. 5957E-08 

4 . 4000E00 

3. 4599E-08 

4 . 3000E00 

9. 0391E-08 

4 . 2000E00 

2 . 3216E-07 

4 . 1000E00 

4 . 4053E-05 

4 . 0000E00 

2 . 5596E-04 

3 . 9000E00 

1 . 2381E-04 

3 . 8000E00 

2 . 6967E-04 

3 . 7000E00 

6. 5175E-04 

3 . 6000E00 

5. 9104E-04 

3 . 5000E00 

3. 7149E-04 

3 . 4000E00 

8. 9525E-04 

3 . 3000E00 

5 . 2663E-04 

3 . 2000E00 

1 . 0270E-03 

3 . 1000E00 

1 . 12  93E-03 

3 . 0000E00 

1 . 1971E-03 

2 . 9000E00 

1 . 1819E-03 

2 . 8000E00 

9. 7390E-04 

2 . 7000E00 

6. 6506E-03 

2 . 6000E00 

6. 0077E-03 

2 . 5000E00 

1 . 72  66E-03 

2 . 4000E00 

7 . 3536E-03 

2 . 3000E00 

1 . 0536E-02 

2 . 2000E00 

6. 7957E-03 

2 . 1000E00 

1 . 0070E-02 

2 . 0000E00 

5. 9813E-03 

1 . 9000E00 

7 . 1253E-03 

1 . 8000E00 

9. 9440E-03 

1 . 7000E00 

1 . 1682E-02 

1 . 6000E00 

1 . 3600E-02 

1 . 5000E00 

4 . 8977E-02 

1 . 4000E00 

2 . 1962E-02 
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1 .3000E00 
1 .2000E00 
1 . 1000E00 
1 . OOOOEOO 
9 . OOOOE-Ol 
8 . OOOOE-Ol 
7 .OOOOE-Ol 
6. OOOOE-Ol 
5. OOOOE-Ol 
4 .OOOOE-Ol 
3. OOOOE-Ol 
2 .OOOOE-Ol 
1 .OOOOE-Ol 
0. OOOOEOO 


1. 6226E-02 
2 . 2041E-02 
4 . 22  62E-02 
5 . 3605E-02 
1 . 1609E-01 
5 . 357  9E-02 
5 . 3704E-02 
6 . 2564E-02 
8 . 4903E-02 
1 . 2800E-01 
5 . 1 98  6E-02 
8 . 4751E-02 
5 . 2  631E-02 
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9 0 % / 1 0 %  U-235/U-238  at  1  Hour 


Photon  Source  Definition  -  generated  from  ORIGEN  F71  file 
Filename:  today. f71 

case  106.  3.600E03  seconds  0.000E00 

total  strength:  2.2393E19  gammas/second 
average  energy:  0.8220924456853023 

discrete  lines:  0.0000E00  gammas/second 
0.00%  of  energy 
average  line  energy:  0.0 

multigroup  bins:  2.2393E19  gammas/second 
100.00%  of  energy 

average  group  energy:  0.8220924456853023 

not  counted:  5.5186E16  gammas/second 

1.10%  of  the  original  energy 


Multigroup  bins  (in  MeV)  and  their  probabilities 


groups:  35 

fraction:  1 

3 . 5000E00 

3. 6600E-04 

3 . 4000E00 

9. 5663E-04 

3 . 3000E00 

5. 7590E-04 

3 . 2000E00 

6. 1879E-04 

3 . 1000E00 

1 . 11 95E-03 

3 . 0000E00 

1 . 2714E-03 

2 . 9000E00 

1 . 2060E-03 

2 . 8000E00 

9 . 2172E-04 

2 . 7000E00 

7 . 1114E-03 

2 . 6000E00 

8.17  65E-03 

2 . 5000E00 

1 . 9391E-03 

2 . 4000E00 

1 . 0378E-02 

2 . 3000E00 

1 . 2126E-02 

2 . 2000E00 

8 . 8413E-03 

2 . 1000E00 

1 . 0719E-02 

2 . 0000E00 

5. 4976E-03 

1 . 9000E00 

9. 5165E-03 

1 . 8000E00 

1 . 0375E-02 

1 . 7000E00 

8. 9979E-03 

1 . 6000E00 

1 . 3411E-02 

1 . 5000E00 

5. 4511E-02 

1 . 4000E00 

3. 0795E-02 

1 .3000E00 

2 . 0595E-02 

1 .2000E00 

2 . 1935E-02 

1 . 1000E00 

5. 1278E-02 

1 . 0000E00 

5. 5379E-02 

9 . 0000E-01 

1 . 1571E-01 

8 . 0000E-01 

5 . 2353E-02 

7 . 0000E-01 

5 . 1749E-02 

6 . 0000E-01 

5 . 6240E-02 

5 . 0000E-01 

9 . 3161E-02 

4 . 0000E-01 

8 . 4548E-02 

3 . 0000E-01 

5 . 6112E-02 

2 . 0000E-01 

8 . 6115E-02 

1 . 0000E-01 

5 . 5395E-02 

0 . 0000E00 
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Excerpt  of  Isotope  Contributions  of  90%/10%  U-235/U-238  at  72  Hours 


Photon  Source  Definition  -  generated  from  ORIGEN  F71  file 
Filename:  today. f 71 

case  125.  2.592E05  seconds  0.000E00 

total  strength:  1.3112E17  gammas /second 
average  energy:  0.42279692832610744 

discrete  lines:  1.3112E17  gammas/second 
100.00%  of  energy 

average  line  energy:  0.42279692832610744 

multigroup  bins:  0.0000E00  gammas/second 

0.00%  of  energy 
average  group  energy:  0.0 

not  counted:  0.0000E00  gammas/second 

0.00%  of  the  original  energy 


Discrete  lines  (in  MeV) 

1 . 5900E00  1 . 2  662E-31 
1.5914E00  8 . 0662E-08 
1.5918E00  1 . 3275E-26 
1.5925E00  6 . 7550E-22 
1.5927E00  4 . 7  924E-1 9 
1 . 5930E00  1.8146E-14 
1.5931E00  3 . 1148E-05 
1 . 5932E00  2 . 5 694E-0 6 
1.5938E00  2 . 8235E-14 
1.5945E00  2 . 4707E-31 
1.5947E00  4 . 6602E-26 
1.5953E00  6 . 6337E-20 
1 . 5960E00  1.7296E-14 
1 . 5961E00  3 . 0011E-23 
1 . 5962E00  2.5568E-02 
1 . 5964E00  3 . 3 600E-37 
1 . 5965E00  3.2226E-14 
1 . 5966E00  8 . 9195E-28 
1 . 5967E00  1.4203E-07 
1 . 5967E00  7 . 3566E-28 
1 . 5968E00  1 . 1 032E-24 
1.5969E00  2.0042E-34 
1 . 5973E00  3 . 4591E-14 
1.5976E00  1 . 5713E-21 
1 . 5985E00  1 . 8  830E-2  6 
1 . 5985E00  2 . 5184E-08 
1.5993E00  5 . 8 683E-2 6 
1.5996E00  1 . 1 983E-30 
1.6000E00  1.7528E-14 


and  their  probabilities 

Dy-155 

Sn-125 

Tm-1 66 

Sr-83 

Ag-llOm 

Eu-154 

1-132 

Sb-128 

Pa-234 

Dy-155 

Bi-214 

Ba-139 

Cd-117 

Pr-140 

La-140 

Er-161 

Eu-154 

Pr-139 

Ga_  7  2 

Tm-1 66 

As-72 

Eu-158 

Cd-117 

Sr-83 

As -7  1 

Sb-129 

Bi-214 

Dy-155 

Sn-127 
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Excerpt  of  Isotope  Contributions  of  90%/10%  U-235/U-238  at  1  Month 


Photon  Source  Definition  -  generated  from  ORIGEN  F71  file 
Filename:  today. f 71 

case  137.  2.630E06  seconds  0.000E00 

total  strength:  7.2408E15  gammas /second 
average  energy:  0.629771890136701 

discrete  lines:  7.2408E15  gammas/second 
100.00%  of  energy 

average  line  energy:  0.629771890136701 

multigroup  bins:  0.0000E00  gammas/second 

0.00%  of  energy 
average  group  energy:  0.0 

not  counted:  0.0000E00  gammas/second 

0.00%  of  the  original  energy 


Discrete  lines  (in  MeV)  and  their  probabilities 

1 . 5841E00  5 . 4  932E-24  Lu-172 
1 . 5854E00  2 . 6928E-13  Pa-234 
1 . 5882E00  3 . 1849E-21  Ac-228 
1.5899E00  1.0086E-14  1-133 
1 . 5914E00  2 . 0312E-07  Sn-125 
1.5925E00  9 . 3175E-27  Sr-83 
1.5927E00  8 . 0421E-18  Ag-llOm 
1 . 5930E00  3 . 2 662E-13  Eu-154 
1.5931E00  1 . 4  908E-0  6  1-132 
1 . 5932E00  4 . 6439E-27  Sb-128 
1 . 5938E00  6 . 4 627E-13  Pa-234 
1.5947E00  4 . 0811E-22  Bi-214 
1 . 5961E00  1 . 9242E-24  Pr-140 
1.5962E00  1 . 5807E-01  La-140 
1 . 5965E00  5 . 8 004E-13  Eu-154 
1 . 5967E00  1.5376E-10  Ga-72 
1.5968E00  3 . 8346E-29  As-72 
1.5976E00  2 . 1673E-26  Sr-83 
1 . 5985E00  3 . 1341E-28  As-71 
1 . 5993E00  5.1390E-22  Bi-214 
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Appendix  D 

Weighted  Average  Build  Up  Factor  Data  for  90/10  U-235  at  One  Hour 


Incident 

Energy 

(MeV) 

Post- 

Scatter 

Energy 

(MeV) 

%  Air 

Contribution 

Above 

Detector 

%  Air 

Contribution 

Below 

Detector 

%  Contribution 

from  Air 

Scatter 

(Isotropic) 

%  Contribution 

from  Ground 
Scatter  (Isotropic) 

%  Total 

Contribution  2nd 
Flight  (Isotropic) 

GNAT 

Weight 

BUF 

SNC  (J-mA2 
/kg/sec/Bq) 

0.020 

0.019 

61.86% 

38.14% 

35.26% 

1.58% 

36.84% 

0.46 

0.11 

1.463 

3.6186E-16 

0.030 

0.028 

69.98% 

30.02% 

26.55% 

3.69% 

30.24% 

0.44 

0.36 

1.375 

1.6419E-16 

0.040 

0.037 

72.97% 

27.03% 

22.81% 

6.05% 

28.86% 

0.43 

0.27 

1.341 

9.952E-17 

0.050 

0.046 

74.35% 

25.65% 

20.67% 

7.98% 

28.65% 

0.41 

0.04 

1.319 

7.515E-17 

0.060 

0.054 

75.16% 

24.84% 

18.90% 

9.13% 

28.03% 

0.39 

0.01 

1.298 

6.6817E-17 

0.070 

0.062 

75.70% 

24.30% 

17.67% 

9.69% 

27.36% 

0.38 

0.00 

1.281 

6.6291E-17 

0.080 

0.070 

76.13% 

23.87% 

16.61% 

10.04% 

26.65% 

0.37 

0.19 

1.265 

6.9925E-17 

0.090 

0.078 

76.48% 

23.52% 

15.89% 

10.20% 

26.09% 

0.35 

0.02 

1.252 

7.5976E-17 

0.100 

0.086 

76.79% 

23.21% 

15.30% 

10.28% 

25.58% 

0.34 

0.01 

1.241 

8.3765E-17 

0.200 

0.153 

78.86% 

21.14% 

12.20% 

9.35% 

21.55% 

0.25 

0.10 

1.177 

1.8911E-16 

0.300 

0.211 

80.16% 

19.84% 

10.74% 

8.47% 

19.21% 

0.20 

0.06 

1.147 

3.0781E-16 

0.400 

0.260 

81.10% 

18.90% 

9.62% 

7.70% 

17.32% 

0.16 

0.10 

1.126 

4.1981E-16 

0.500 

0.305 

81.84% 

18.16% 

8.77% 

7.10% 

15.87% 

0.13 

0.11 

1.111 

5.3146E-16 

0.600 

0.344 

82.55% 

17.45% 

8.10% 

6.53% 

14.63% 

0.11 

0.06 

1.100 

6.4279E-16 

0.700 

0.380 

82.96% 

17.04% 

7.48% 

6.14% 

13.62% 

0.10 

0.06 

1.090 

7.4302E-16 

0.800 

0.413 

83.39% 

16.61% 

6.95% 

5.73% 

12.67% 

0.09 

0.06 

1.083 

8.5302E-16 

0.900 

0.444 

83.78% 

16.22% 

6.51% 

5.39% 

11.90% 

0.08 

0.13 

1.076 

9.3949E-16 

1.000 

0.472 

84.21% 

15.79% 

6.18% 

5.09% 

11.26% 

0.07 

0.06 

1.072 

1.0327E-15 

1.100 

0.499 

84.41% 

15.59% 

5.77% 

4.81% 

10.58% 

0.06 

0.06 

1.066 

1.1209E-15 

1.200 

0.524 

84.68% 

15.32% 

5.46% 

4.57% 

10.03% 

0.06 

0.02 

1.062 

1.2075E-15 

1.300 

0.547 

84.94% 

15.06% 

5.17% 

4.34% 

9.51% 

0.05 

0.02 

1.059 

1.2913E-15 

1.400 

0.569 

85.18% 

14.82% 

4.92% 

4.13% 

9.05% 

0.05 

0.03 

1.055 

1.3736E-15 

1.500 

0.590 

85.39% 

14.61% 

4.70% 

3.94% 

8.64% 

0.05 

0.06 

1.052 

1.4527E-15 

1.600 

0.610 

85.59% 

14.41% 

4.48% 

3.77% 

8.25% 

0.04 

0.02 

1.050 

1.5289E-15 

1.700 

0.629 

85.77% 

14.23% 

4.29% 

3.61% 

7.90% 

0.04 

0.01 

1.047 

1.6027E-15 

1.800 

0.647 

85.94% 

14.06% 

4.12% 

3.46% 

7.57% 

0.04 

0.01 

1.045 

1.6744E-15 

1.900 

0.665 

86.10% 

13.90% 

3.96% 

3.32% 

7.28% 

0.03 

0.01 

1.043 

1.7439E-15 

2.000 

0.682 

86.25% 

13.75% 

3.81% 

3.20% 

7.01% 

0.03 

0.01 

1.041 

1.8117E-15 

2.500 

0.757 

86.90% 

13.11% 

3.21% 

2.66% 

5.87% 

0.03 

0.00 

1.034 

2.1275E-15 

3.000 

0.822 

87.39% 

12.61% 

2.77% 

2.28% 

5.05% 

0.02 

0.00 

1.029 

2.4073E-15 

3.500 

0.878 

87.80% 

12.20% 

2.44% 

1.97% 

4.41% 

0.02 

0.00 

1.026 

2.6601E-15 
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Weight  Values  From  the  ORIGEN  Fallout  Analysis  Tool  for  90/10  U-235  at  12  hours 


12  Hour  Case  Sample  Calculations 
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